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Abstract 

This lecture note starts with a pedagogical introduction to the theoretical 
background and properties of perfect actions, gives some details on topology 
and instanton solutions and ends with a discussion on the recent developments 
concerning chiral symmetry. 



1 Work supported in part by Schwcizcrischer Nationalfonds and by Iberdrola, Ciencia y Tec- 
nologia, Espafia. 



1 Introduction 



This lecture note is a summary of the authors contribution at different schools in 
1997: 'NATO AS1 - Confinement, Duality and Non-perturbative Aspects of QCD' 
in Cambridge, 'Advanced School on Non-Perturbative Quantum Field Physics' at 
Peniscola and 'Yukawa International Seminar on Non-Perturbative QCD' in Kyoto. 
It is extended by discussing in some detail recent developments on chiral symmetry. 

This lecture note is intended to be largely self-contained. The general aspects 
of lattice regularization and Wilson's renormalization group and in particular 
those points which are closely related to our subject will be explicitly introduced 
and discussed. The reader might even consider our topic as a pretext to learn 
about these powerful theoretical and practical tools in quantum field theories. I 
hope, however, that the unexpected and amazing properties of the perfect actions 
will also catch the readers' fantasy. 

The very definition of a quantum field theory requires a regularization. As it is 
well known, all the regularizations break some of the symmetries of the underlying 
classical theory. They introduce a new scale (the cut-off) which is reflected in the 
final predictions even after the regularization is removed: naive dimensional analysis 
breaks down, anomalous dimensions are created, dynamical mass generation might 
occur, etc. On the other hand, beyond these exciting physical effects, the new 
scale creates unwanted cut-off artifacts everywhere which have to be eliminated in 
a theoretically and practically difficult limiting process. On the tree level (in the 
classical theory) the new introduced scale creates artifacts only, nothing else. If 
the classical theory is scale invariant (like the Yang-Mills theory or massless QCD) 
the cut-off artifacts destroy the most interesting and relevant classical properties of 
the model like the existence of scale invariant classical solutions or fermionic zero 
modes and index theorem. 

Regularizations often break other symmetries also. The lattice, the only known 
non-perturbative regularization, breaks Euclidean rotation symmetry (Lorentz sym- 
metry) which is less of a headache, and also chiral symmetry which creates a difficult 
tuning and renormalization problem. Chiral symmetry awakes deep theoretical is- 
sues (anomalies, no-go theorems) culminating in problems concerning chiral gauge 
theories. 

Perfect actions || offer a bold approach to these problems. We shall call a lattice 
regularized local action classically perfect if its classical predictions (independently 
whether the lattice is fine or coarse, whether the resolution is good or bad) agree 
with those of the continuum. The quantum perfect action does the same for all the 
physical questions in the quantum theory. That such actions exist might seem to 
be surprising. Their existence is closely related to renormalization group theory. 
As we shall see, they have beautiful properties which one would not expect a lattice 
action can have. 

The classically perfect action has scale invariant instanton solutions, has no topo- 
logical artifacts, satisfies the index theorem concerning the fermionic zero modes on 
the lattice, preserves all the important physical consequences of chiral symmetry in 
the quantum theory and is expected to reduce the cut-off effects significantly even 
in quantum simulations. 

In this note we shall mainly discuss the theoretical properties of the classi- 
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cally perfect actions. These actions are the fixed points of renormalization group 
transformations in asymptotically free theories and are determined by classical field 
theory equations. They are not abstract theoretical constructions only. Actually, 
most of the effort during the last four years was directed towards constructing and 
parameterizing them explicitly and testing their performance in numerical simula- 
tions ||]. We shall not discuss these developments here at all beyond the remark 
that the fixed point action gives very good results in quantum simulations includ- 
ing questions related to topology. The main practical question is whether one can 
find a parameterization which approximates the fixed point action sufficiently well 
without making the simulation prohibitively expensive. 



2 Lattice regular izat ion and the continuum limit 

In this Section we summarize some of the basic notions concerning lattice regular- 
ization, continuum limit and the relation between Quantum Field Theories (QFTs) 
and critical phenomena in classical statistical mechanics. Readers, who are familiar 
with these basics are advised to move directly to Section 2.4 on locality. 



2.1 QFTs defined by path integrals in Euclidean space 

In these lectures we shall consider QFTs in their Euclidean formulation and use path 
integrals. Using path integrals to describe quantum theories not only provides an 
appealing physical picture for the quantum evolution as a sum over classical paths 
(in quantum mechanics) or sum over classical configurations (in QFT), but gives 
a powerful framework for analytic manipulations and opens the way for numerical 
calculations also. In addition, it helps to understand the deep relation between 
QFTs and critical phenomena in classical statistical mechanics [Q, |5| . 

The formal relation between QFTs and classical statistical physics is easy to see. 
Consider an n-point function in a scalar QFT in Minkowski space 

(0\T(fa{)fa 2 )...$(x n j)\0}, (1) 

where 4>{x) is the field operator. In the path integral language this expression is 
given by 



D(()<f>(xi)(f>(x 2 )...4>(x n )exp(i / dtdx C(<f>, d^)), (2) 



1 

~Z 

where Z = J Dcj) exp(i J dtdx C) provides for the correct normalization, x = (t, x) = 
(t, X\, . . . , Xd s ), and the Lagrangian is defined as 

1 ds 1 

i=i 

The relation we are looking for becomes clear if we go to Euclidean space by rotating 
time to imaginary time 

t = x Q —> -ix d (4) 

It can be shown B that in every order of perturbation theory the analytic con- 
tinuation in eq. (0) (Wick rotation) is possible due to the specific pole structure 
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of the Feynman propagators. It is an assumption that this remains true beyond 
perturbation theory. Using 

i J dtdx -> J d d x, (5) 

and 

d 1 

C-Mink -» -(^ -^d^d^cj) + V (</>)) = -C-Euch (6) 

we get 

Z = y" D<£exp(- /" d d x£ £ „ c/ (0,cV)), (7) 



(<p(xi) . . . (f>(x n )} = — I Dcj)(t)(xi)...(f>(x n )exp(- I d d xC E ud)- (8) 

These equations can be interpreted as the partition function and the correlation 
function of a system in classical statistical mechanics. The quantum theory of fields 
in d s space dimensions is transformed into classical statistical mechanics of fields in 
d = d s + 1 Euclidean dimensions. The Euclidean action plays the role of E/kgT, 
where E is the classical energy of the d dimensional configuration. 

Actually, the relation between QFTs and classical statistical systems is deeper 
and more specific than the observation above. We shall return to this question in 
Sect.2.3 



2.2 Lattice regularization 



Some of the mathematical operations which enter the definition of a QFT require a 
careful limiting procedure. In field theory, the variables are associated with space- 
time points. As can be seen from the form of the Lagrangians of different classical 
FT's (sec, for example eq. (^)), these variables have some kind of self- interaction, 
whereas the elementary interaction between different degrees of freedom is over in- 
finitesimal distances as expressed by derivatives. Already in the classical theory, the 
definition of a derivative requires the temporary introduction of a finite increment 
(of the argument of the function considered) which disappears at the end by some 
limiting procedure 

f {x) = um f±±^LJM . (9) 

a^Q a 

Defining the derivative, which has the important role of providing communication 
between different field variables, is a much less trivial problem in QFTs. Like 
in eq. (^), the very definition of a QFT requires the temporary introduction of a 
defining framework called regularization, which disappears from the theory by a 
limiting process. The way to introduce and remove the regularization is a highly 
non-trivial problem. 

We shall use lattice regularization in these lectures. A hypercubic lattice is 
introduced in the d = d s + 1 dimensional Euclidean space. Scalar and fermion field 
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variables live on the lattice points, vector fields on the connecting links. Derivatives 
are replaced by some kind of finite difference operation. 

It will be convenient to work with dimensionless quantities. In the following, 
all the fields, momenta, masses, couplings and other possible parameters will be 
defined dimensionless by absorbing appropriate powers of the lattice unit a in their 
definition. For any quantity Qdf with dimension s (measured in mass units) we 
define 

Qdf -> Q = a s Qdf ■ (10) 

For a scalar field 4>{x) we write 

4>{x)^(j> n , (li) 

where n is a d-dimensional vector of integers. The simple choice for the finite 
difference operation 

d^ix) -> V M 0„ = (j> n+ p, - 4> n (12) 
leads to the following lattice regularized form of the classical Euclidean Lagrangian 

A= [ d d x£^Y.J2\ V ^nV^n+Y. V (M- ( 13 ) 

n /i n 



Having a finite lattice unit a the lattice defines a short distance regularization. 
Going over to momentum space, in the Fourier integral exp(ipn) enters which is 
periodic under — > p M + 2tt. This constrains the momentum to the Brillouin zone 
— t > Pn < I"- Therefore, the lattice provides a momentum cut-off also. 

The lattice regularization has several advantages over other conventional reg- 
ularizations used in perturbation theory. In general, in problems which require 
numerical analysis (e.g. integrals, differential equations,...), the standard procedure 
is to introduce meshes. Also the path integral can be defined in a natural way by 
introducing meshes. The path integral regularized this way becomes a set of inte- 
grals and is ready for numerical procedures in case of non-perturbative problems. 
Actually, lattice regularization is the only known non-perturbative regularization. 
It is a very natural regularization also for theories with gauge invariance or with 
constrained variables (like the non- linear tr-model). 

Undoubtedly, lattice regularization has some disadvantages also. It breaks cer- 
tain symmetries, as every regularization does. But the lattice has special problems 
(although not unsurmountable) with chiral symmetry and tough problems with chi- 
ral gauge theories. The discretization of space-time breaks Lorentz symmetry (0(d) 
symmetry in Euclidean space), but, as we shall see later (Sect. 8.1), the remaining 
cubic symmetry is sufficient to have it restored automatically as the regularization 
is removed (at least in asymptotically free theories). The lack of Poincare symmetry 
creates problems, however, in theories with supersymmctry. 

Although the lattice is used, in general, to investigate non-perturbative prob- 
lems, for certain questions (renormalization of operators, finding improved actions 
or comparing lattice results with those obtained in other regularizations) it is un- 
avoidable to use lattice perturbation theory. Unfortunately, perturbation theory on 
the lattice is technically somewhat cumbersome. 

Some aspects of lattice regularization will be discussed in Sect. 5. For further 
reading we refer to textbooks . 
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2.3 The continuum limit of QFTs versus critical phenomena 
in statistical mechanics 

Assume, we can solve our lattice regularized QFT at some set of (bare) couplings 
and masses in the Lagrangian. The predicted spectrum will contain dimensionless 
numbers, since all our quantities were made dimensionless by absorbing appropriate 
powers of the lattice unit a (Sect. 2.2). The mass of an excitation is given by a 
number M, the corresponding dimensionfull mass Mdf, [Mdf] = g and correlation 
length tdf, [fa] = cm are given by 

M dS = M\ t* = j^ = ±* = t*- (14) 

The correlation length defined above is the length scale over which the massive 
particle can propagate with a significant amplitude. For a generic value of the 
couplings and mass parameters in the Lagrangian the predicted M will be an 0(1) 
number, the mass of the particle is of the order of the cut-off 1/a, the correlation 
length is O(a). This corresponds to a very poor resolution, the situation is far from 
the continuum limit. It requires a careful tuning of a certain number of parameters 
in the Lagrangian to make the predicted M's in the low lying spectrum to become 
very small numbers which correspond to physical length scales much larger than 
the lattice unit a. In this limit the lattice becomes very fine and the presence of the 
rcgularization will not distort the physical results (like mass ratios) anymore. This is 
the process of removing the lattice regularization which leads to the continuum limit. 
Theories in which this can be achieved by tuning a finite number of parameters (and 
by fixing the normalization of fields appropriately) are called renormalizablc. This 
criterion is identical to that introduced in perturbation theory. 

In the continuum limit the correlation length of physical excitations is much 
larger than the lattice unit a. But how large is it in Fermi? This can not be pre- 
dicted, the absolute scale should come from observations. If the theory is expected 
to describe the physics of hadrons, for example, then the lowest lying particle with 
spin=l/2, electric chargc=l, baryon numbcr=l, strangcncss=0,... should be identi- 
fied with the proton having a mass of 940 MeV. This fixes the physical length scale 
and relative to that unit the lattice constant goes to zero in the continuum limit. 

In Sect. 2.2 we observed a formal relation between QFTs and classical statistical 
systems. We see now that the continuum limit of a QFT corresponds to a statistical 
system with a correlation length which is much larger than the lattice unit a. This 
is the case for critical statistical systems. In critical solid state problems the lattice 
unit is fixed (typically 0(A)) and the correlation length becomes macroscopical. 
This is a matter of choosing the absolute scale, however, and it does not reduce the 
extent of the deep analogy between the two subjects. 



2.4 Locality 

A recurrent issue in our discussion will be the locality of the action density. The 
laws of classical physics (from Newton through Maxwell to Einstein) are expressed 
in terms of differential equations (rather than, say, integro-diffcrcntial equations). 
They correspond to actions, where the interaction between variables in different 
space-time points extends over infinitesimal distances as expressed by derivatives. 
This continues to be true in quantum systems also. No symmetry principles would 
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prevent us to add a term to the Lagrangian in eq. like 

C{x) = ... + J d d y f(y)<f>(x + V -)<P(x ~ |) , (15) 

where f(y) is an even function of y and has an extension A (for example f(y) ~ 
exp(— y 2 /2\ 2 )), where A is finite when measured in Fermi, or it might decay only 
as a power in y 2 . Nevertheless, such models found their place neither theoretically, 
nor experimentally. 

In a QFT with non-local interactions renormalizability (in the sense discussed 
in Sect. 2.3) will be lost. In addition, a vital and beautiful property of local QFTs 
(shared by critical statistical systems), called universality will be lost also. Univer- 
sality means that the physical predictions become independent of the microscopic 
details of the Lagrangian, i.e. they are not sensitive to the detailed form of the 
interaction at the cut-off scale. Since our knowledge on interactions much above 
the scales of present experiments is very limited, universality is vital to preserve the 
predictive power of a QFT. 

On the lattice we shall call an action density local, if it has an extension of 0(a). 
Typically, it is an exponential function of the distance |n — n'\ between the variables 
at n and n': ~ exp(— j\n — n'\), 7 = 0(1). In the continuum limit (a — > 0), the 
extension of the action density measured in physical units goes to zero. On the 
other hand, an action density which has a finite width in physical units is non-local. 
No such non-local actions will be considered here in the following. 

A Lagrangian which has nearest-neighbour interactions only or interactions 
which arc identically zero beyond a few lattice units is certainly local, but no phys- 
ical principle requires to have this extreme case. If this were necessary to assure 
the above mentioned nice properties of a QFT (or of a critical statistical system) , 
then no experiment on ferromagnets close to the Curie temperature would observe 
universality. Really, the elementary interaction between the magnetic moments in 
the crystal decays rapidly, but certainly does not become identically zero beyond a 
few lattice spacings. 

In the following we shall call an action where the couplings become identically 
zero beyond a certain 0(a) range ultralocal. We shall also use the expression 'short 
ranged' for local actions where the exponential decay goes with a large 7. 



3 Renormalization group 

A QFT is defined over a large span of scales from low (relative to the cut-off) physical 
scales up to the cut-off which goes to infinity in the continuum limit. Although 
field variables associated with very high scales do influence the physical predictions 
through a complicated cascade process, no physical question involves them directly. 
Their presence and indirect influence makes it difficult to establish an intuitive 
connection between the form of the interaction and the final expected predictions. 
The presence of a large number of degrees of freedom makes the problem technically 
difficult also. It is, therefore a natural idea to integrate them out in the path integral. 
This process, which reduces the number of degrees of freedom, taking into account 
their effect on the remaining variables exactly, is called a renormalization group 
transformation @, § (RGT). 
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In this Section we introduce some of the basic notions related to RG theory and 
set the notations. Since these lectures are mainly concerned with asymptotically 
free (AF) QFTs, our discussion will be biased by this goal. 



3.1 Renormalization group transformation 

For simplicity consider a scalar field theory regularized on the lattice in d Euclidean 
dimensions as discussed in Section 2.2. The RGT averages out the short distance 
fluctuations, i.e. the fluctuations over distances O(a). The lattice is indexed by 
n, the field associated with the point n is denoted by </>„. We introduce a blocked 
lattice with lattice unit a' = 2a whose points are labeled by ub and the associated 
block field is denoted by Xn B ■ The block variable Xn B 1S an average of the original 
4> fields in the neighbourhood of the point Ub- 

Xn B = b ^2 ^{2n B - n)(j) n , (16) 

n 

where u) specifies the weights of the averaging whose normalization is chosen as 

^w(2n B -n) = l. (17) 

n 

The significance of the scale factor b in eq. ([!(]) will become clear later. A trivial 
choice for lo might be to take uj{2ub — n) = 2~ d if n is in the hypercube whose 
center is indexed by ns and zero otherwise, fig. 1. 




Figure 1: A simple example for a block transformation in d=2. 
We integrate out now the original <p n variables keeping the block averages fixed: 
exp(-^'(x)) =11 / #"11^^ - bJ2^( 2n B - n)<j) n ) exp(-A(4>)) . 

n ns n (18) 

The new action A'(x) describes the interaction between the block variables. The 
partition function remains unchanged 

Z ' = II / d ^ B exp(-.A'(x)) = J] / exp(-.A(0)) = Z . (19) 

The long-distance behaviour^ of Green's functions, and so the spectrum and other 
low energy properties of the system is expected to remain unchanged as well. On 

2 'Long-distance' is meant here and elsewhere in this work as 'many lattice units', and it might 
mean short, or long physical distances. 
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the other hand, the lattice unit is increased by a factor of 2: a! = 2a. Since the 
dimensions are carried by the lattice unit (Sect. 2.2), the dimensionless correlation 
length, which is measured in the actual lattice distance, is reduced by a factor of 2: 

£ - | ■ (20) 

By eliminating the fluctuations on the shortest scales we reduced the number of 
variables by a factor of 2 d , while taking their effect into account by changing the 
action A — > A' appropriately. Although the lattice unit is increased by a factor 
of 2, and so the resolution became worse, the long-distance behaviour remained 
unchanged. In particular, no new cut-off effects are generated in the predictions 
for physical quantities. Iterating this RGT, the goal formulated at the beginning of 
Sect. 3 is achieved. 

We close this subsection by generalizing the RGT of eq. (|l^) slightly || (a step 
which will be very useful later): 

exp(-A'(x)) = 

JV{k)Y[ / d<f> n exp[-A(<j)) - Kj2(Xn B -6^w(2n s -n)0„) 2 ], (21) 

n ub u 

where k is a free parameter. For k — + oo eq. (21) goes over to eq. (|l8|), for finite k 
the block average is allowed to fluctuate around x (slightly, if k is not too small). 
Therefore, the parameter k specifies the stiffness of the RG averaging. N is a 
normalization factor (which is field independent, trivial in this case) introduced to 
keep eq. ( |l9| ) valid. 



3.2 Constraints on the block transformation from 
symmetries 

It is largely arbitrary how the block averages in a RGT are constructed, but the 
procedure should conform with the intuitive goal of a RGT: it should lead to the 
elimination of the short distance fluctuations. Identifying, for example, the block 
variable with one of the original fine variables in the block ('decimation') is not 
really an averaging, and although it is a legal transformation of the path integral, 
it will not lead to useful results in general. 

Beyond this general requirement, symmetries put further constraints on the form 
of the block transformation. Denote the fields on the fine and coarse lattice by <f> 
and x? respectively and write the RGT in the form 

exp(-A'(x)) = I D <t> II /(XWccfe;0)exp(-^)). (22) 

blocks 

Here (f>, x can be scalar, vector or fermion fields and / defines how the coarse field 
in the block Xbiock is constructed from the (neighbouring) <fi fields. 

Assume, there is a symmetry transformation <j> —> R(f> under which the action 
and the measure are invariant: A{R(f>) = A((j>), D{R(j>) = D<f>. The blocked action 
A'(x) wm inherit this symmetry, A'(Rx) — ■A!(x)i if the averaging function satisfies 

f(Rxbiock; 4>) = fixbiock^R- 1 ^ , (23) 
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as it is seen easily by changing integration variables in eq. ( p2| ) . 

In case of gauge symmetries we shall require somewhat more: we shall not only 
require that the effect of a gauge transformation on the coarse lattice is equivalent to 
a gauge transformation on the fine lattice, eq. (^3|), but also the other way around. 
This constraint implies that gauge equivalent fine configurations contribute equally 
to the path integral in eq. (p2|). 



3.3 A basic assumption of the RG theory 

Consider a step of the RGT (eq. ( |l8| ) or eq. (21)) for the case where the parameters of 
A(<p) were chosen so that the original system was close to the continuum (or, in the 
language of statistical physics, close to criticality). As eq. (21) shows, performing a 
RGT on a system with action A(<fi) leads to a path integral where A is replaced by: 

A(0) + K J2(Xn B -bYs^ns-n)^) 2 . (24) 

riB n 

For this path integral the field \ enters as an external field. Even if A((f>) was 
carefully tuned to criticality, the system in eq. (|2~i| ) is not expected to be critical. 
Actually, this is a basic assumption of the RG theory: the path integral entering 
a step of the RGT defines a non-critical problem with short-range interactions 
only. This is what one expects intuitively: the \ fields constrain the (/>-averages 
in their neighbourhood and disrupt the long-range fluctuations. This assumption 
leads to the important conclusion that the new action A'(x) will be local. Really, 
if the r.h.s. of eq. (21) describes a system with short-range fluctuations only, the 
generated interaction between distant x fields is expected to be negligibly small. 

Another consequence of this basic assumption is that the path integral in eq. ( |l8| ) 
or eq. (21) is a technically much simpler problem than the path integral of the 
original system. Saying differently: to perform a RG step is not an easy problem 
in general, but it is much simpler than to solve the original (near) critical theory. 
This is what makes the RG theory a practically useful idea. 



3.4 Fixed Point (FP) and the behaviour 
in the vicinity of a FP 

In general, the transformed action A'(x) in eq. (21) will contain all kinds of interac- 
tions, even if the original action A(4>) had a simple form. It is useful to introduce a 
sufficiently general interaction space to describe the actions generated by the RGT, 
and write 

A(<f>) = J2 K M4>), (25) 

a 

where 9 a ((j)), a — 1,2,... denotes the different interaction terms like (for scalar 
fields) 

01 = V P<A« > ^ = > • • • > h = Yl ^V^nV^ n , . . . 

n,fi n n,fi (26) 



10 



and K a are the corresponding dimensionless couplings. The transformed action 
A'(x) is also expanded in terms of these interaction terms 

A'( X ) = J2 K 'MX)- (27) 

a 

The RGT induces a motion in the coupling constant sp&ce:{K a } — > {K' a }. Under 
repeated RGTs a coupling constant flow is generated 

{*r } (28) 

while the (dimensionless) correlation length is reduced with every step (eq. (EOT)) 



(29) 



It might happen that certain points in the coupling constant space are repro- 
duced by the RGT 

{K*} -> {K*} . (30) 

A point with this property is called a fixed point (FP) of the RGT and the corre- 
sponding action 

a fp w = J2k*M) (31) 

a 

is the FP action. FP actions play an important role in QFTs in general and will 
play an important role in our later discussion of the perfect actions. 

Eq. (^8|) and eq. (^) imply that at the FP, £ = oo or £ = 0. We will be interested 
in FPs with £ = oo. The set of points in the coupling constant space where £ = oo, 
forms a hypersurface, which is called the critical surface. As eq. (|2^) shows, an 
RGT drives the point {K a } away from the critical surface, except when {K a } is on 
the critical surface. 

Let us consider now the behaviour of the flow under the RGT in the vicinity of 
a FP 0. Take a point {K a } with AK a = K a - K* small. Under the RGT we 
have 

AK a -» AK' a =K' a -K* a . (32) 
Expanding AK' a = AK' a ({K}) around {K*} 



AK' a = AK' a ({K*}) + J2 j^-AK' a ({K}) lK=K ,AK p + 0({AKf) , 

dK P (33) 



and observing that the first term on the r.h.s. is zero, we get the following linearized 
RG equation 

AK> a = J2T a pAKp, (34) 





with 



T a0 = ^AK' a ({K}) {K=K , . (35) 
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Let us denote the eigenvectors and eigenvalues of the matrix T by h a and A a , 
respectively]^] 

Y,T a ph a f} = X a h a a ,a = l,2,... . (36) 

P 

The eigenvectors h a define the eigenoperators 

ft°(0=5>°0«(0)- ( 3? ) 

a 

These operators define a new basis in terms of which the action can be expanded 

A(0) = A FP (<b) + Y,c a h a W, (38) 

a 

where c a are the corresponding expansion coefficients (couplings) which are small 
if A is close to A FP '. Repeated application of the RGT gives 

A{cj>) - A FP {<f) + A°c /i°(^) -> A FP {<f) + ^(A Q )V/i a (0) - . . . , 

(39) 

i.e., the coupling c a goes over to (A a )™c a after n RGT steps. For |A a | > 1 (|A a | < 1) 
the coupling \c a \ is increasing (decreasing) under repeated RG steps. The corre- 
sponding interaction h a is called relevant (irrelevant) . For A a = 1 the fate of the 
operator (called marginal) is decided by the higher order corrections suppressed in 
eq. @. 



3.5 RG flows of a Yang-Mills theory 

The following discussion applies also for other AF theories, like the O(N) non-linear 
a- model or the CP n model in d = 2. 

The classical action of an SU (N) Yang-Mills theory has the form 

PA cont = 4tJ x F ^ x)F ^ W ' a = " • • ' N2 ~ 1 ' (40) 

where (3 = 2N/g 2 and F^ v is the colour field strength tensor. Any lattice represen- 
tation of this action should go over to the form in eq. (|4^) for slowly changing weak 
fields. Even if we use a gauge symmetric discretization (as we always do), there 
is a large freedom in writing down a local action with this property on the lattice. 
As we discussed in Sect. 3.4, we have to introduce an infinite dimensional coupling 
constant space if we want to follow the change of the action under RGTs. We write 

[3Ai att = PAi aU (K u K 2 ,...), (41) 

where we indicated the coupling constant dependence of the action only. Under a 
RGT we have: {/?, Ki,K 2 , ...} -> {p , K[,K' 2 , ...}. 

The expected flow diagram is sketched in fig. 2. In the j3 = oo hyperplanc, 
which is the critical surface, there is a FP with coordinates K*,K%, ■■■■ We shall 

3 Actually, the real matrix T is not symmetric in general. For a discussion on related conse- 
quences, see the first reference in [nj 
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Figure 2: Flow of the couplings under RG transformations in a Yang-Mills gauge 
theory. 

suppress the index 'latt' and introduce the notation A(K^, K2, ■■■) = A FP . All the 
eigenoperators which lie in the (3 = 00 hypersurface are irrelevant. There is one 
marginal direction which is pointing out of this surface. Actually, this direction is 
marginal in the linear approximation only, it becomes weakly relevant by the higher 
order corrections to eq. (p3|). 

Consider the point {(3, K\, K 2 , ...}, (3 very large. Under RGTs this point runs 
rapidly towards the FP (the largest irrelevant eigenvalue is 1/4), then slowly moves 
away from it along the weakly relevant direction tracing a trajectory, called the 
renormalized trajectory. 

Consider the interaction corresponding to the FP, A FP , multiply it by (3: (3A FP 
and allow (3 to move away from 00. This is a straight line (the dashed line in fig. 2.) 
which is not a RG flow, but defines an action for every value of (3. We shall call 
this action the FP action. 

We postpone the discussion of the arguments supporting the structure of the 
flow diagram in fig. 2. until Sect. 8. We mention here only that theories with this 
flow diagram can be renormalized by tuning a single coupling constant and the 
resulting renormalized theory is universal (independent of the value of the other 
couplings in the (bare) action). The continuum renormalized theory is obtained as 
g — > (/? — > 00) in which limit the correlation length goes to infinity, the lattice unit 
a — > 0, the resolution becomes infinitely good and the cut-off artifacts disappear 
from the predictions. 
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4 Perfect actions in AF theories 



We shall call a lattice regularized local action classically perfect if its classical pre- 
dictions (independently whether the lattice is fine or coarse, whether the resolution 
is good or bad) agree with those of the continuum. The quantum perfect action does 
the same for all the physical questions in the quantum theory. That such actions 
exist might seem to be surprising. As we shall see, they do exist and have beautiful 
properties which one would not expect a lattice action can have [§, ||[ |ll|, [j^, [l3|, [l4| . 

We shall argue now that the actions defined by the points of the renormalized 
trajectory (RT) in fig. 2 define quantum perfect actions. It means that by taking 
an arbitrary point ((3, K\, K2, ■••) G i?T, where g 2 — 2N//3 need not be small, 
the corresponding action A(/3, Ki^K^,, ...) will produce quantum results for physical 
questions (mass ratios, for example) which are exactly the same as in the continuum 
limit. This is surprising, since in the point (/3, K\, K2, ■■■) the correlation length is 
not large, the lattice unit a is not small, the resolution is not fine. Nevertheless, 
the statement is true. 

The argument goes as follows^ . At any given /3, the point of the RT is connected 
to the infinitesimal neighbourhood of the FP by (infinitely many steps of) exact 
RG transformations. Since each step increases the lattice unit by a factor of 2, any 
distance at the given f3 (even 1 lattice unit) corresponds to a long distance close 
to the FP. The infinitesimal neighbourhood of the FP is in the continuum limit, 
there are no cut-off effects there at long distances. On the other hand, for all the 
questions which can be formulated in terms of the degrees of freedom after the 
transformation we get the same answer as before the transformation. Thus, there 
are no lattice artifacts at the given j3 on the RT at any distances. 

We shall also argue that the FP action, as defined in Sect. 3.5 (fig. 2) is a clas- 
sically perfect action. We are not able to demonstrate this statement so compactly 
as that for the quantum perfect action before. In the following sections we shall go 
over the different classical properties of continuum Yang-Mills, QCD and other AF 
theories and show one by one that they are reproduced by the FP action indepen- 
dently of the lattice resolution. An additional, unexpected property of the QCD 
FP action is, as we shall discuss in some detail, that it reproduces the physical 
consequences of chiral symmetry in the quantum theory exactly. 

The detailed form of the FP action and the RT depend on the form and param- 
eters of the block transformation. The theoretical properties of the FP action and 
RT are, however, independent of these details. In this sense, any RGT which has 
a local FP is suitable. On the other hand, different local FP actions (belonging to 
different block transformations) have different extensions and this has an important 
practical significance when it comes to parameterization and simulation. In order 
to find short ranged actions (Sect. 2.4), different block transformations were tested 
and optimized || . 
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5 Lattice regularization and explicit RGTs in dif- 
ferent AF theories 



It is time now to come down from generalities to concrete examples. We shall 
discuss how the lattice regularization is set up and give explicit RGTs in different 
AF theories. First, the 0(A) non-linear cr-model and the SU(N) Yang-Mills theory 
will be considered. Then, as a preparation for QCD, the special difficulties of lattice 
regularized Dirac fermions will be illustrated by treating the case of free spin-1/2 
fermions. We close this section with QCD. 

In all these theories we write down the 'standard' lattice action which is the 
simplest discretization of the continuum action satisfying the basic requirements. 
These are the actions which - due to their technical simplicity - were used in most of 
the numerical simulations. These actions are simple, but they represent poorly most 
of the properties even of the classical field theory if the resolution is not sufficiently 
fine. As we discussed in Sect. 4, in the framework of the RG theory actions can be 
defined, which have beautiful properties. For later use we define explicit RGTs in 
all the theories considered in this section. 



5.1 O(N) non-linear cr-model in d=2 

The classical action in the continuum has the form 



where S is an A-component vector satisfying the constraint S 2 (ir) = 1, (3 = 1/g. 
Writing 



and expanding in the fluctuations ir(x) (it has A — 1 components), a systematic 
perturbation theory can be set up. The theory is AF in the coupling g for A > 3. 
It is believed (and this is strongly supported by different theoretical and numerical 
results) that a non-zero mass is generated dynamically, the spectrum contains a 
massive O(N) multiplet. The a- model is discussed in different modern textbooks 
on QFT and statistical systems || O, |lG| . 

Looking for classical solutions with finite action, the spin vectors S should be 
parallel at large distances. For A = 3 this leads to an S2 — > 5*2 mapping between 
the coordinate and group spaces. The 0(3) model has a non-trivial topology, there 
exist scale invariant instanton solutions with an action A = Air\Q\, where Q is the 
topological charge of the solution. 

The 0(3) model shows many analogies [l6| with the Yang-Mills theory in 
d = 4. It offers the possibility to test non-perturbative ideas and numerical methods. 

On the lattice, the scalar field S lives on the lattice points. The simplest real- 
ization of the action has the form ('standard action') 




(42) 




(43) 




(44) 
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where V M is defined in eq. ( p"2| ) and = 1. On smooth configurations A s t goes 
over to Acont ~ a property we shall require for any lattice action. 

Consider RGTs with a scale factor of 2, forming averages out of the four fine 
spins in a block. The block spin will be denoted by R„ B , R^ B = 1- The RGT has 
the form 

cxp(-/3U'(R)) = J DSexp[-[3(A{S) +T CT (R,S))] , (45) 
where T a defines the averaging process and we used the notation 

J DS = IlJ dS >AS 2 - 1) . (46) 



We shall give two explicit examples for the averaging which were studied in the 
context of perfect actions g, |l9). Both RGTs are O(N) invariant, satisfying eq. ( f23[) 
for global O(N) transformations. In the first case: 

T CTl (R,S) = ^(- Kff R„ B £„ B +JV(|£„ B |)) , (47) 

TIB 

where E nB is the sum of the four S spins in the block ub 

£„ B = £ S„ , (48) 

n£n B 

and Af assures the correct normalization, eq. (|l9|). We get 

M(\H nB \) = ~hiY N (pK a \Il nB \) , (49) 

where Y^(z) is related to the modified Bessel function (some of its properties are 
summarized in appendix F in poj , for example), specifically Y^(z) ~ sin(z)/z. For 
n a — > cxd, the block transformation goes over to a (^-function constraining the average 
to lie parallel to the coarse spin 

v 

R„ B , (k„ = oo) . (50) 



It will be useful, however, to keep k g finite and optimize its value to get a short 
ranged FP action. 



Modifying slightly the averaging function T 

S 



T ff2 = 2 KCT ^(R„ B -^) 2 , (51) 



tlB I 



we can make the technically cumbersome normalizing function Af constant. Really, 
the integral 

J dRS(K 2 - 1) exp(4/3K CT eR) , (52) 

with e = S/|S|, depends on the length of e only, therefore it is an S-independent 
constant. 
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5.2 SU(N) Yang-Mills theory in d=4 



The classical action in the continuum has the form 

f3A™(A») = wj x Tr ( F n" F ^) > (53) 

where = 2N/g 2 , = F£ u T a , a = 1, TV 2 - 1 and T a are the generators of the 
colour group SU(N), Tr{T a T b ) = l/25 ab . 

The classical theory has scale invariant instanton solutions [fl8| with an action A = 
&t: 2 \Q\. 

The quantum theory is AF in the coupling g. It is believed (and this is supported 
by many, but non-rigorous results) that non-perturbative effects create a confining 
potential between static sources in the fundamental colour representation (static 
quarks) and the spectrum contains massive colour singlet excitations (glucballs) 
only. 

On the lattice, the vector gauge field lives on the links of the hypercubic lattice 
]2l] , 0. The convenient variables, in terms of which gauge invariance is easily kept 
on the lattice are not the vector potentials themselves, but their exponentiated 
forms. Consider a static quark and antiquark in the points x and y, respectively 
(in the continuum). The form 

^(i)Pexp(t f A„{z)dz^) tP{y) , (54) 

where the integral runs along some path between x and y and V denotes path 
ordering, is gauge invariant. The gauge part in eq. (|54|), which leads the colour 
flux between the sources, is an element of the SU(N) group rather than that of the 
algebra. Going to the lattice and placing x, y in the endpoints of the link (n, n + fi), 
we get 

i>nUn{n)%l>n+p,, (55) 

where U^n) is related to the lattice vector potential as 

U^n) = expiiaA^n)) . (56) 

The variables of the lattice regularized Yang-Mills theory are the link matrices 
Ufj,(n). The link matrix associated with the same link but directed oppositely: 
(n + ft, n) is f/T(n). Under a gauge transformation the link matrix transforms as 

Up(n) -» V(n)U^n)V^(n + fi) , (57) 

where V(n) € SU(N) is the gauge transformation in the point n. Consequently, 
the trace of the product of link matrices along a closed path on the lattice is gauge 
invariant. 

The simplest lattice realization of the Yang-Mills action has the form ('Wilson 
action') 

PA W (U) = jjJ2( N - ReTr U p) > ( 58 ) 

plaq 
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where the sum is over the plaquettes and U p is the product of the four directed 
link matrices around the plaquette p. For smooth gauge fields eq. ( p8| ) goes over to 
eq. (|53|), as it should. 



We shall consider RGTs with a scale factor of 2. The blocked link variable 
V^ub) G SU(N), which lives on the coarse lattice with lattice unit a' = 2a, is 
coupled to a local average of the original link variables. The blocking kernel T(V, U) 
enters the RGT as in the c-model 



eiq>(-f3'A'(V))= I DU exp[-p(A(U) +T g (V,U))] , (59) 

where 

/ DU = Iij dU M- ( 6 °) 



The kernel T g is taken in the form 



T g <y,U) = £ [-^ReTr(V^n B )Ql{n B )) +N{Q»{n B ))} . (61) 



N 

n B ,M 



The N x N complex matrix Q^ns) is an average of the fine link variables in 
the neighbourhood of the coarse link l B — (n B ,/j,). N{Q n(n B )) is fixed by the 
normalization condition, eq. (|19|): 



exp(7V(Q)) = J dWexp(^-ReTr(WQ' i )) : W G SU(N) . (62) 

Let us give a simple example for the construction of the average link matrix Q 
HID (%• 3): 



QnM = (1 - 6c) [^(n^n + A) + 

c [Uv{n)U^{n + P)U fl (n + A + 0)Ul{n + 2fi)+ 

Ul{n - v)U„(n - v)Uf,{n + fi - v)U v {n + 2fi - v)] , (63) 

where c, the relative weight of the staples versus the central link, is a tunable 
parameter (see, last paragraph of Sect. 4). 



' i 3 • ' = (l-6c) B > O > O + c* o + ... 

Figure 3: A simple example for a gauge invariant blocking 

The averaging process in eq. (63) is not a very good one: there are many link 
variables on the fine lattice which do not contribute to any of the Q's. Experience 
shows that poor averaging might lead to a FP which is local, but not sufficiently 
short ranged or might not even have a FP. It is not difficult, however, to generalize 
the block transformation above and to make it physically more satisfactory p3[ . 
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5.3 Free fermion fields in d=4. Chiral symmetry, doubling 
and no-go theorems 



The action of spin 1/2 Dirac fermions in the continuum has the form 

^cont(-0» = / -0(^)(7m 5 m + m)ip(x) , (64) 

J x 

where {7^,7^} = 25^ and 75 = 71727374- 

On the lattice, the fermion field lives on the lattice points. The most general 
action which is invariant under cubic rotation, reflection, permutation of the coor- 
dinate axes and discrete translations, reads 

A{$, ip) = ^ ^Hri - n')ip n > , (65) 

n,n' 

where 

h(n — n 1 ) = p^{n — n')^^ + X(n — n') . (66) 
The local functions p M and A have the following reflection properties 

ft»(-r) = -p„(r), A(-r)=A(r), (67) 

and so h satisfies 

h) = l 5 h l5 . (68) 
A naive discretization of eq. ( |64]) leads to 

Anaiveii), VO = \ ^ {^nl^n+jl ~ ^n+£7^n) + ^ m ^>n^n , (69) 



2 



corresponding to 



Pu(r) = ^(Sr+p,,o - 5r-p,,o), A(r) = m8 rfi , (70) 
or, in Fourier space 

Pfj,(k) =isin(fc M ), A(fc)=m. (71) 

The momentum space propagator has the form 

— z y^,, sin(fc u )7„ + m 

D(k) = s V M7 M . (72) 

E M sm (fc M ) + m 2 

The poles of the denominator determine the spectrum. Since sin(fc^) = at = 
and fc M = 7r, the propagator in eq. ( [72]) has 16 poles defining 16 fermion species 
with the correct relativistic dispersion relation E = —ik^ — (k 2 + m 2 ) 1 / 2 in the 
continuum limit. This is the species doubling problem specific to fermions on the 
lattice. 

For ttl — 0, the action in eq. ( |64| ) is chiral invariant, i.e. it remains unchanged 
under the transformation 

tp(x) exp(ia75) ip(x) , ip(x) — » ip(x) exp(i«75) . (73) 
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This symmetry requires on the lattice 



{h,Tt>} = 0, (74) 

where we denoted the anticommutator by {,}. The naive action e g. (|69| ) satisfies 
cq. ([74]) for m = 0. In general, eq. (|7|) require A = in cq. (|66|) and leads 
to a propagator Pfiln) _1 - F° r small fc, p^(k) ~ ifc^ since the lattice action 
should reproduce the continuum action in this limit. In addition, p^(fc) is a periodic 
function (Sect. 2.2). It follows then that —ip^k) > for small and positive and 
< for k^ close to, but less than 2ir. If p^(k) is a continuous function, it should 
have a zero in between. This intuitive argument implies a close relation between 
the chiral symmetry of the action and the doubling problem: if the action is chiral 
symmetric, we have species doubling The no-go theorem says: if the lattice 

action has the correct continuum limit, i.e. h(k) = ij^k^ for small momenta and 
there are no other zeros in the Brillouin zone (no doublers), then one of the following 
two conditions should be violated: a/ h(k) is an analytic periodic function of 
b/ {^75} = 0. Clearly, we do not want to violate condition a/ because this would 
imply loosing locality. 



So, we have to violate the chiral symmetry of the action. It is well known 
in perturbation theory in models without chiral symmetry (scalar QFT) or when 
non-chiral invariant regularization is used (QED, for example, with Pauli-Villars 
regularization) that radiative corrections induce an additive mass renormalization 
and so, a tuning problem. This will be the generic situation in our case also. Simi- 
larly, the axial current will not be conserved, the soft pion theorems will be violated 
and there will occur mixing between operators in different nominal chiral represen- 
tations. However, as we discuss later, the FP action offers an elegant solution to all 
these problems. 

A simple way to kill the doublers is to add to the action an appropriate non- 
chiral invariant term which disappears in the continuum limit. Wilson suggested 
l25[] ('Wilson fermion action') to add to eq. ( |69|) 

- ^ (tyni'n ~ i>n4>n+fi ~ ^n+fi^n) (75) 

with r^Oa free parameter. This term contributes to \{k) 

\{k) = m + 2r^sin 2 (^). (76) 

The extra term gives an 0(1/ a) mass to the doublers, but it is suppressed around 
the k ~ pole. This procedure can be generalized to interactive theories also and 
this is the method which is used in most of the QCD calculations. 

Consider now RGTs on a d-dimensional hypercubic lattice with a scale factor of 
2. Write the RGT in the form 

exp(-^'(x,x)) = [ D$Dil>exp[-(A$,4)+T f (x,xrf,il>))] (77) 



with 



J DtPDtP = II / d ^n#n , (78) 
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where ip,ip are Grassmann variables. Consider blocking kernels of the form 
Tf(x, X', ip) = K f ^2 (Xn B ~ b f^2 i> n u(2n B - n)) 

x (xn B -bfY, Lo{2n B -n)i> n ) . (79) 

n 

We shall assume that the averaging function uj is diagonal in Dirac space and real. 
Observe that Tf is not chiral invariant. 

We shall consider two different block transformations j^, |2?| . The first is 
just the blocking introduced for a scalar field in Sect. 3.1 and also for the cr-model 
in Sect. 5.1: lo(2ub — n) = 2~ d if n is in the hypercube whose center in indexed by 
tib and zero otherwise. 

Keeping in mind QCD, the blocking in the second example has a form which is 
easy to make gauge invariant in the presence of gauge interactions. Place the block 
fields in the even points of the fine lattice (n is even if all its coordinates are even). 
A simple choice for the averaging function u) in this case is 

w(0,0,... ,0) = c , w(±l,0,... ,0) = ... = u;(0,0,... ,±l) = d, 

... > w(±l,±l,...,±l) = c a . (80) 

The parameters cq, c\, ...,Cd satisfy the normalization condition, eq. (fl7|) ; 

c + 2da + 2d(d - 1) c 2 + . . . + 2 d c d = 1 . (81) 

Observe that the fine field ip n contributes to several block averages, in general. We 
shall call a blocking with this property 'overlapping'. The block transformations 
in eqs. ( (47}|5l| ) and the first example above are 'non-overlapping', while those in 
eqs. (63,80) are overlapping. 

There is a natural symmetric choice for the a parameters in eq. (80): 

c t = 2- d -\ i = 0,l,...,d. (82) 

In this case, the sum of contributions of ip n to all the block fields x is independent 
of n, the averaging is 'flat'. 



5.4 QCD 

The action is the sum of the Yang- Mills action and the action of quarks in interaction 
with the colour gauge field 

Jx (83) 

The quark fields carry Dirac, colour and flavour indices, the gauge field does not 
know about flavour, the mass (matrix) does not know about colour. For N = 3, 
this theory is believed to be the theory of the hadrons. It should explain the 
masses, widths and other static properties of hadrons, scattering events, all kinds of 
hadronic matrix elements, the behaviour of hadronic matter at finite temperature 
and density, spontaneous chiral symmetry breaking, the rj' problem, etc. Most of 
these problems are non-perturbative. 
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The simplest lattice regularized QCD action, which we shall call the 'Wilson 
action' is the sum of the gauge action eq. ( |5S| ) and the fermionic action, eq. ( |69| ) 
plus eq. (frq), made gauge invariant. Take r = 1, which is the preferred value (for 
technical reasons) in simulations 

A% CD tl>, U) = A £(JV - i?eTr f/ p ) + 

2^ft(> " l)E^(r# n +/J - &,+a(7/* + l)Ut(n)ih) + ^ ^„(m + 4)^„ . (84) 

n,fi n 



We write the RGT in the form 

exp[-(/3X(F) + 4(x,X,F))] = 

D$Dil>DUexp[-0(Ag{U) + T g (V, U)) + {Aj®, i>, U) + T f ( x , x; & V>, U))] ■ 

(85) 

The gauge kernel T g is defined in eq. (|6l|), while the fcrmion kernel is the gauge 
invariant version of eq. (79) 

T f(X, X, i>, i>, U) = Kf ^ (Xn B ~ b f^2 $ n uH U )n,n B ) 

TLB n 

x (Xn B - b f^ 

riB ,n i>n). (86) 

n 

In eq. (86), the averaging function lu is made gauge invariant by connecting rig with 
the point n by a product of J7 matrices. This can be done without problems if tie 
overlaps with one of the fine lattice points as in the second example in Sect. 5.3, 
eq. (80). Beyond hypercubic symmetry the choice of these paths is part of the free- 
dom we have in defining an RGT. A simple choice is to take the shortest connecting 
paths. For U ->• 1, oj(U) 

riB ,n — > lu(2i~ib — n) defined in eq. (79), which satisfies the 
normalization condition, eq. (|l7j). This condition gives a meaning to bf in eq. (86). 
As we shall see (Sect. 8), in a free field theory bf is determined by the engineering 
dimension of the field. This should be the case also in the classical theory even in 
the presence of interaction, since the corrections to bf are related to the anomalous 
dimension of the field generated by quantum fluctuations |8|, || . We shall return 
to this point later. 



6 The saddle point equation for the fixed point 
action 



For free field theories with a quadratic blocking kernel, like those in eq. (79) or in 
eq. (21), the RGT is reduced to Gaussian integrals. The integrals can be performed 
exactly and the problem can be investigated in every detail. In general, however, 
the path integral of a RGT is a highly non-trivial problem. 

A basic observation to proceed is that in AF theories the FP lies at (3 = oo and 
in this limit the path integral can be calculated in the saddle point approximation 
g]. This leads to an equation in classical field theory which determines the FP 
action. We shall derive this equation first in the non-linear a-model and then in 
QCD. 
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6.1 The FP of the O(N) non-linear a- model 



For a given fixed coarse field configuration R we can perform the path integral 

y"^Sexp[-/3(^(S)+r ff (R,S))] (87) 

in the (3 — > oo limit using the saddle point approximation leading to the classical 
equation 

A! (R) = min [A(S) + T2° (R, S)] , (88) 

where is the j3 = oo limit of the blocking kernel. The FP of the transformation 
is determined by the equation 

.4 PP (R) = min[^ p (S) +T»(R,S)] . (89) 

Let us write down this equation explicitly for the block transformation defined in 
eqs. ©|§): 

.4 PP (R) = min[.4 pp (S) - Kct ^(R„ b S„ b - |£ nB Q] , (90) 

TLB 

where we have used lnljv(z) = z(l + 0(\nz/z)) for large z and S„ B is defined in 
eq. ©. 

Observe that eq. ((8?]) is reduced to the saddle point equation, eq. (^8[) , for any 
configuration R. If the configuration R is strongly fluctuating then the minimizing 
configuration S m i n will not be smooth either. In general, eqs. ( |S8| , [39| ) and their 
solutions have nothing to do with perturbation theory. A starting condition for 
the FP equation, eq. (^9|), is that A FP goes over to the classical action ^4 co „ t (R), 
eq. (|4^), for very smooth R configurations. 



6.2 The FP of QCD 

In the f3 — > oo limit (V,\,X fixed), the Boltzmann factor on the r.h.s. of eq. (85) 
is dominated by the gauge part. The integral over the U field is saturated by the 
minimizing configuration U m in — U m in(V) leading to the equation 

A' g (V) = mm[A g (U) + T~(V, U)] . (91) 

The gauge part of the FP is determined by 

A FP (V) = mm[A FP (U)+T~(V,U)] . (92) 

Using the explicit form of the blocking kernel T g in eq. ( |6l| ) we obtain 

A FP (V) = mm (a FP (U) - ^ £ [RzTr(V»MQl(n B )) ~ /(<?>*))]] , 

m /(93) 



where 



f(Q) = max\ReTr(WQ^)] , W G SU(N) . (94) 
w 
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Eq. ( p4[ ) follows from eq. ( p2| ) in the f3 — > oo limit. 
The remaining fermionic integral has the form 

exp(-A' f (x,X,V)) = / D^I?V'exp[-(^ / (^,V',^mm)+7>(x I X;^V' ) C / mm))], 
J (95) 

where [/ m j n (7) is the minimizing configuration from eq. (pl|). If the action Af is 
quadratic in the fermion fields, then the blocked action A'j is quadratic also, since 
the kernel Tf is quadratic (eq. (86)) and the integral in eq. ( |95| ) is Gaussian. Write 
the fermion action in the following general form 

Af{i>^,U)=Y,i>nK, nl {U)i)' n , (96) 

n,n' 

where h(U) n ^ n i carries Dirac, colour and flavour indices which are not indicated 
explicitly. Using eq. (86) and performing the Gaussian integral in eq. (^5|) we get 
the following recursion relation for h 

h (V)n B n' B = ^5n B n' B (97) 

-Kffrf ^2^{U mi n)n B n (h{U min ) + Kfbf {U min )uJ {U min )) ^ w(Z7 roin )t, n> . 
nn' 

In case h~ x is defined (no zero modes, see Sect. 10) eq. ( |97| ) is equivalent to the 
somewhat simpler equation p6[ p8[ 

h '( V )nl,n' B = ^*»B,n' B +^Sw(tUi)ns,«ft(^n)^i l /W(Knia)^, ni • 

The FP satisfies the equation 

hFP (VV B w = — <W B + b 2 f u(U mia ) nBtn h FP (U min )-^ n ,uj(U mia )l, „, . 

' B K f n,n> ' B (99) 

Eqs. ( p3| , p9| ) determine the FP action of QCD. Gaussian integrals of c-number fields 
are equivalent to minimization. This can be generalized to Gaussian Grassmann 
integrals to express eq. (§9) as 



(100) 



where ip st ,ip s t make the r.h.s. stationary, 5/8ip [...]=0, S/Sip [...]=0. In eq. ([10 1 
the fields x, x, t/>, ?/> are c-number fields. 



7 The FP action for weak fields 



We have claimed in Sect. 4 that the FP action is classically perfect: it reproduces 
all the essential physical properties of the continuum classical action on the lattice 
even if the resolution of the lattice is poor. The proof of this statement follows 
directly from the FP equations (eq. (|9^) for the cr-model or eqs. ([)3],[)9|) for QCD) 
without solving them explicitly. Before presenting these arguments let us study the 
FP equations in a limit where analytic results can be obtained. 
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Consider first the FP equation in the O(N) non-linear cr-model, eq. (|90[). Take 
a configuration where the spins fluctuate around the first axis in O(N) space 



where x ns has (N— 1) components. Assume, these fluctuations are 'weak', Xn B I *C 
1. In this case, the minimizing field S m i„(R) will also fluctuate around the first 
axis, eq. (p3| ), and these fluctuations will also be weak, \ir nB \ <C 1. For weak 
fields eq. ( |90| ) can be expanded in powers of the fluctuations and one can study the 
solution order by order. 

There is a difference between the 'weak field' introduced here and the 'smooth 
field' used repeatedly earlier. A smooth field changes slowly in configuration space 
and is dominated by small k values in Fourier space. In a weak field, the size 
of deviations from the classical vacuum configuration is small, but the field might 
change rapidly in configuration space, so higher momentum components might be 
important in Fourier space. 

In leading order of the weak field expansion, the FP equation is quadratic in the 
fluctuations. Using translation symmetry we can write 

A FP (K) = ±Y1 P(r B )Xn B X nB+ r B + 0(X 4 ) , (102) 

ns-XB 

where the unknown couplings, to be determined by the FP equation, are denoted 
by p{tb)- Expanding also the blocking kernel in eq. (|90|), we get at the quadratic 
level 



g P( r B)Xn B Xn B +r B = (103) 

UB,TB 

Tx^k ^ P( r )*mrn+r + 2 n a ^ (Xn B ~ \ 7r "^] ' 



Using the trivial relation between Gaussian integrals and minimization we observe 
that eq. (Co3) is just the FP equation of a free scalar theory: 



CeX P(^ H P( r B)Xn B Xn B +r B ) = (104) 
n B 

J[ / d7r n exp[--^p(r)7r„7r„ +r -2«; .^(x nB - ^ ^ 7r„) 2 ] . 

n n, r n B n£n B 

The same steps can be used for the Yang-Mills or QCD FP actions. In quadratic 
order the Yang-Mills FP equation, eq. (|93]), goes over to the FP problem of a free 
(Abelian) gauge theory, while the fermion part of the QCD FP action, eq. (p^), 
becomes equivalent to the FP equation of free Dirac fermions. In the next section 
we give a general solution of the free field FP problem. 



8 The free field FP problem 



As we discussed in Sect. 7, in leading order of the weak field expansion the FP 
equations of AF theories are reduced to that of free fields. We shall study the free 
field problem in this section. 
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In all our previous considerations we applied lattice regularization. For free 
field theories without gauge symmetry we might also use a simple momentum space 
cut-off, which leads quickly to some results we already referred to in the previous 
sections]!], ^, 30 . We present these results first, and then turn to the problem of 



FP actions on the lattice. 



8.1 Scalar field with momentum space cut-off 

We work in Fourier space and constrain all the momenta to be smaller than A cut . As 
always, we use dimensionless quantities, the dimensions are carried by the cut-off. 

We write the quadratic action in the general form 

A{<j>) = \( <K-q)p(q)<l>(q), (105) 

and perform a RGT by integrating out the fields in the momentum range 1/2 < 
\q\ < 1. We write 

<M<z) = (106) 

where </>o(<?) and 4>\{q) are non-zero only for \q\ < 1/2 and 1/2 < |g| < 1, respec- 
tively. In the RGT we have to integrate out 4>\{q), i.e. perform the path integral 



exp(-A(<f> )) = II / ex P 
l<kl<i 



(A((h + <k))- (107) 



Using the quadratic action in eq. (105) we get 



Afa + fa) = ~ [ <f>(-q)p{q)<f>(q) + \ [ H-q)p(q)Hq) ■ 

1 J \q\<k ZJ i<\i\<i (108) 



Writing eq. fllOq ) into eq. ( |105| ), the (/^-dependent part can be brought out of the 
integral, while the 0i-integral gives a field independent, irrelevant constant. By 
rescaling and relabeling the 4>o field as 

W) = C"Vo(?) 9' = 2<?, 0<|q'|<1, (109) 

we get 



AM=C 2 \l <f>'(-2q)p(q)0'(2q) 



2 



? 2 - d \ [ 4/(-q)p(q/2)<l/(q)=A'W). (110) 
2 J\q\<i 



The role of the rescaling factor £ in eq. (109) is the same as that of b in eqs. (Mjlq) 



whose significance we shall see in a moment. As eq. (105) and eq. (110) show, under 
the RGT p transforms as 

p(q) - C 2 2- d p(q/2) . (Ill) 
Expand p in powers of momenta 

p(q) =a + a 2 ^2 + a4 + «6 ^m) 3 • • • ■ 



(112) 
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For ao small, the pole of the propagator l/p(q) is close to q — and ao — m 2 , 
where m is the mass of the free particle. We shall require that the action density 
before and after the transformation goes over into the form of classical field theory 
for small q- values, i.e. we fix the coefficient of q^q^i in p to be 1. This condition 
fixes C 

C 2 2- d Eyy-E^-C = 2^. (iia) 



Eq. (Ill) gives then the transformation rules for the other couplings 

a — > 2 2 ao , a 4 — > 2 -2 a 4 , ccq — > 2~ 3 ag , . . . . (H4) 



Eq. (114) shows that the operator <fi 2 (the coupling ao or m 2 ) is relevant with an 
eigenvalue Ao = 4, all the other operators (couplings) which contain four or more 
derivatives are irrelevant. The largest irrelevant eigenvalue is 1/4. 

Replacing the region of allowed momenta of a cut-off sphere by a cut-off hyper- 
cube, the regularization will violate rotation symmetry. Only hypercubic symmetry 



remains, like on the lattice. In this case we have to allow in eq. (112) terms like 
A ^fiilp-Qv-Qp-Q^)' which is not 0{d) invariant, but hypercubic symmetric. The cou- 
pling /3 4 is also irrelevant, /3 4 — ► l/4/?4 under the RGT. That means that rotational 
symmetry is restored as the RGT is iterated, as we go towards the infrared. 

If we start with ao = 0, we are and remain on the critical surface and run rapidly 
towards the FP which, in this case, has the simple form 



A FP (^) = ll c)>(-q)q^(q) , (lir,| 

kl<i 



2 



describing a free massless scalar particle. As we mentioned before, the form of the 
FP depends on the block transformation (typically, it will look more complicated 
on the lattice), but the conclusions concerning the relevant, marginal and irrelevant 
eigenvalues is independent of that. 



Consider now small perturbations around the FP, eq. (115), like A 4 </> 4 , A 6 </> 6 , • 
with A 4 , Ag, • • • small, as discussed in Sect. 3.4 

A(4>) = A FP {^) + ^ f 4>(-q)cP(q) + A 4 / <K<h) • ■ ■ ^)S(qi + . ■ ■ + g 4 ) 



A 6 / <j>{qx) . . . <t>{q G )8{qi + . . . + q t ) + . . . , (116) 

Jqi,— ,96 

where only three of the perturbations are written out explicitly. Using this action A 
in eq. ( |107| ), expanding the exponent in ao, A 4 , \& . . . up to linear order, performing 



the Gaussian o -mtegrals and rescaling and relabeling as in eq. (109), we get 



a -» (a + 6cA 4 + 15c 2 A 6 + ...)C 2 2~ <i , 

A 4 -» (A 4 + 15c 2 A 6 + ...)( 4 2- 3d , (117) 

A 6 -» (A 6 + ...)C 6 2~ M , 

with 




(118) 
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The corresponding T-matrix (eq. (35)) is a triangular matrix, the eigenvalues A are 
just the diagonal elements 



A 1= 2 2 , h<»= / 2 (x), 

J X 

A 2 = 2 4 - d , h& = [ [^(x) + acj) 2 (x)] , (119) 

J X 

A 3 = 2 6 ~ 2d , h<® = f [<jfi( x ) + fc<f{x) + (3 2 <j> 2 {x)] . 



As eq. (119) shows, there is a simple relation between the eigenvalue Xg and the 
dimension dg of the operator 9(x) which has the highest dimension in the corre- 
sponding eigenoperator h: 

Xg = 2 d - de (120) 

In the eigenvector h the highest dimensional operator 9(x) is mixed with lower 
dimensional operators. 

In d = 4, operators with dimension less than, equal, larger than 4 are relevant, 
marginal and irrelevant, respectively. This result is also valid in gauge and fcrmion 
theories. In SU(N) Yang-Mills theory the lowest dimensional gauge invariant op- 
erator is Tt(F^F^), which is marginal. 



8.2 Lattice FP actions in different free field theories 

Consider first a real scalar field. Denote the fields on the coarse and on the fine 
lattice by x an d 4>, respectively. The lattice units are ao and a\ = do/2. Assume, 
we are in the FP, and perform a RGT: 



CeX P(-^ P( r B)Xn B Xn B +r B ) = 

d(/>„exp[--^p(r)0„0„ +r - -K^2(xn B - &^w(2n B - n)<j> n ) 2 ] . (121) 

rt n,T n B n 

Multiply this equation by Xn' B Xri^ an d integrate over the \ fields. Using eq. ( |l9| ) 
we obtain 

(Xn'Xn") = b 2 Y] u{2ri B - n')u)(2n B - ra")(0 n '0n") + ~<*n' «" ■ 

n',n" K (122) 



In Fourier space eq. (|122j) has the form 

p(q B ) 2 i 2 P{2 +7Tl ) 2 K 

l = (h,... ,l d ),lv = 0,1. (123) 



/(">= /" ■■• T T^exp^n)/^). (124) 



The Fourier transformation is defined as 

We demand, as discussed in Sect. 8.1, that any action should go over the classical 
continuum action for smooth fields. That requires that p(<?s) — > q% + 0(q B ) for 
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small qs- (The FP action is on the critical surface, the mass parameter ao in 



eq. (112) is zero.) Considering eq. (123) in the qs — > limit, the I ^ contribution 



on the r.h.s. can be neglected and we obtain 

1 1 , 1 



A -" 2 ^ m2 T^JW- <125) 



The normalization condition, eq. (|17]), gives w(0) = 1 and so 

b = 2^. (126) 



We shall use the relation of eq. (122) recursively, i.e. we connect the propagator 



) on the r.h.s. of eq. (122) with the propagator on the lattice with unit 



o-i = «o/4- Denote the field on this lattice by r\ and the lattice points by nf. 

(Xn' B Xn B ) = & 4 ^( 2t 4 - n')uj{2n - n' f ) 



n'n 



f 

E 

n" ,n'. 



u;(2n B - n")uj(2n" - n'}){r] n ' f Vn';) + 
1 



K 



[S n ' B ,n» + b 2 2>(2n' B - n)u;(2n B - n)] . (127) 



Iterating this equation further, after j steps one arrives at a lattice with lattice unit 
dj = clq/2 j . For j — > oo, the lattice will be infinitely fine and the propagator for 
this lattice can be replaced by the propagator in the continuum. This way, the FP 
propagator (x n ' B Xn B ) will be expressed explicitly in terms of the propagator in the 
continuum and of the product and sum of the averaging function u. The inverse 
propagator gives the FP action. 

If the block transformation is non-overlapping (see the paragraph after eq. (|8l|)), 
the sum in the last term of eq. (127) is ~ l//t<J n ' : „" . This remains true in every 
order of the iteration. Although for overlapping transformations this is not true 
anymore, the ~ 1/k contribution to the propagator remains local. The non-locality 
of the propagator is carried completely by the first term in eq. (127). 

Without going through the detailed derivation, we quote the final results only. 
Let us start with an example: take the simple non-overlapping transformation, 
where uj{2ub — n) = 2~ d if n is in the hypercube whose center is indexed by 
tib and zero otherwise (Sect. 3.1). The FP propagator in Fourier space has the 
form§, HI, H 

1 - V 1 TT ^l + A^ , 4 

p(q) ^(^ + 2^)2 11 ( ^+^)2 + 3k < I ») 

where the summation is over integer vectors l = (h, • • -Jd) and (<7+27rZ) 2 = ^ M (g Al -|- 
27r/ At ) 2 . The ~ 1/k term in eq. ( |l2S| ) is a constant in Fourier space, hence ~ 5 n ,n' m 
configuration space, which is the general result for non-overlapping transformations. 



The second factor under the sum in eq. (128) is regular, it is the Fourier transform of 



the product of the local averaging function to. The pole singularities come entirely 
from the (q + 2irl)~ 2 piece. These poles define the spectrum of the FP theory on 
the lattice. Fix the spatial part of the momentum to be q = p,Pi 6 (— n, it), and 
find the pole singularities in the complex g^-plane. There is an infinite sequence of 
poles on the imaginary axis 

q d = ±i |p + 2tt1| , leZ". (129) 
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Table 1: The couplings of the free field FP action at a distance r = (ro, rj.) for the 
optimal choice of the block transformation with k = 4. Note that in our convention, 
for the standard action the only non-vanishing entry in this list would be Pst(1, 0) — 
-1. 



r 


p(r) 


r 


p{r) 


(1,0) 


-0.61802 


(4,0) 


-2.632 • 10~ B 


(1,1) 


-0.19033 


(4,1) 


7.064- 10~ 7 


(2,0) 


-1.998- 10- 3 


(4,2) 


1.327- 10~ 6 


(2,1) 


-6.793- 10- 4 


(4,3) 


-7.953- 10- 7 


(2,2) 


1.625- 10~ 3 


(4,4) 


6.895- 10" 8 


(3,0) 


-1.173- 10" 4 


(5,0) 


-8.831 • 10~ 8 


(3,1) 


1.942- 10~ 5 


(5,1) 


3.457- 10- 8 


(3,2) 


5.232- 10~ 5 


(5,2) 


3.491 • 10~ 8 


(3,3) 


-1.226- 10" 5 


(5,3) 


-3.349- 10~ 8 






(5,4) 


8.408- 10" 9 






(5,5) 


-1.657- 10" 10 



Denoting p+27rl by k, fc, 6 (— oo, oo), we get the exact relativistic dispersion relation 
of a massless particle 

-iq d = E{\ S ) = \k\. (130) 



Through the Z-summation in eq. (|128j ), the full exact continuum spectrum is repro- 
duced (the momentum k is not constrained to the Brillouin zone), even though we 
are on the lattice, which resolves the large k waves very poorly. This result is to 
be compared with the dispersion relation of the standard nearest-neighbour action 
(cq. ( |l3| ) with V = 0). If k is close to the end of the Brillouin zone, the dispersion 
relation is wrong by a factor of 2 or more. 



Performing the sum in eq. (128) and calculating p(r) by Fourier transformation 
is an easy numerical exercise. The function p(r) and hence the FP action is local, as 
expected: it decays exponentially ~ exp(— 7|r|), as discussed in Sect. 2.4. The size 
of 7 which defines how short ranged the action is, depends on the free parameter k. 
The value k r* 4 is optimal. For d = 2, we have 7 « 3.4 for k = 4. In this case p(r) 
is strongly dominated by the nearest neighbour and diagonal couplings as table 1 
shows |Q . 

Actually, p(r) is a perfect discretization of the continuum Laplacian. The ques- 
tion raised by this observation is, whether the FP idea can be used in solving partial 
differential equations, like the Navier-Stokes equation, numerically fl33| . 

Let us present now the general result. The fields can be scalars, fermions or 
vectors. For the blocking kernel we write 



n B n,0 n,/3 (131) 



where the bar above the fields denotes complex or Dirac conjugation^, the index a 

4 For a real scalar or vector field tj> = <f> and we include a factor of 1/2 both in the action and in 
the kernel, as in eq. (121). 
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can be a Lorentz or Dirac index^. The FP propagator is given by 

D FP {q) = lim £>W(q), (132) 

n — >oo 

where 



2™-l 



= (2-V)« £ ^ ( " ) (^)^ (0) (^)^ t( " ) (^) (133) 



n-l 



3=0 



where 17 and Q are defined as 



Co{2- 1 )0j{2- z )...Cj{2- 1 ), 

2 3 -l 



j(i) 



(g) = (^)i£fl«(^)nt«(^) 



The rescaling factor 6 is given by 



i«=0 



6 = 2 d f 



2-i 



2-' 



(134) 
(135) 

(136) 



where is the dimension of the field <j>. For non-overlapping transformations 
Q (?) is independent of g, in the overlapping case it is an analytic function of q. 
For illustration, consider free Dirac fermions in d = 4 using the overlapping block 
transformation defined by the eqs. (80,^). In this casein 



Il cos2 (y)> 



fjW(g) 

= ^ r n[( 1 + 4 ^) C0S ^) + (2 + 4- J )] 



TT2-* sin2(f } 

2-J 



(137) 



For the initial fermion propagator one might take the Wilson fermion propa- 
gator (Sect. 5.3). Eq. (134) and eq. (137) gives the FP fermion propagator 



D FP (q) = -ia»(q) 7fl +l3(q). 



where 



J* (a) = V gM + 2^ TT sin 4 (f) 



and 



/3(g) = ^[ciJJcos(g A1 ) + c 2 (cos(gi)cos((72)cos(q3) + 

+(cos(g 2 ) cos(g 3 ) cos(g 4 )) + c 3 (cos(gi) cos(<? 2 ) + 
c 4 2jcos(g jU ) +c 5 ] . 



(138) 
(139) 

(140) 

cos(g 3 )cos(g 4 )) + 



We shall consider the block transformation only where u) is diagonal in the Dirac indices. 
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The constants ci , . . . , eg result from different geometric sums and have the value 
ci = 0.592581, c 2 = 1.388406, c 3 = 3.349502, c 4 = 8.336553, c 5 = 102.410571. 
This FP propagator has all the features we expected. The spectrum is given by the 
~ part and it is exact. The ^ 1/k part is analytic in g-space and ultralocal 

(Sect. 2.4) in configuration space. The propagator is not chiral symmetric. The 
chiral symmetry breaking term is ~ 1/k comes entirely from the non-chiral invariant 
block transformation. 



9 Instanton solutions and the topological charge 
on the lattice 

In this section we shall consider the 0(3) non- linear a- model in d = 2. The argu- 
ments and results can be immediately generalized to Yang-Mills theories in d = 4. 

As we discussed in Sect. 5.1, the configurations of the 0(3) model fall into 
different topological sectors characterized by an integer number Q, the topological 
charge of the configuration. While configurations from the same topological sector 
can be continuously deformed into each other, this is not possible for configurations 
with a different charge Q. The charge Q is the number of times the internal sphere 
is wrapped as the coordinate sphere is traversed. It may be defined as 

Q = ^J d 2 xe^S{d^S x 0„S) , (141) 

and it is related to the action by the inequality 

A>Air\Q\. (142) 

If, for a given configuration the equality is satisfied, the configuration minimizes the 
action for the given topological charge Q and is therefore a solution of the equations 



of motion. In order to show the inequality in eq. (142) one writes down the identity 



J lJ 2J (143) 

which follows easily from vector product identities and from dnS 2 = 0. The l.h.s. 



of eq. (143) is non- negative, while on the r.h.s. a combination of the action and of 



the topological charge enters leading to eq. (1143) 



9.1 Instanton solutions of the FP action 



Replacing the continuous Euclidean space by a discrete set of points, the notion 
of a 'continuous deformation of the configuration' is lost and the definition of the 
topological charge bec ome s problematic. A possibility is ('field theoretical defini- 
tion') to discretize eq. (141) in some simple way taking care of the surviving discrete 
symmetries . The nice geometrical properties of the topological charge are lost 
this way, it is not an integer valued functional of the fields any more and it suffers 
from a multiplicative renormalization on the lattice [ |35"| . Another possibility is the 
'geometrical definition' p& which assigns an integer charge to each configuration. 
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However, this charge receives, in general, unphysical contributions from topological 
artifacts as discussed below. 

An additional problem is that the lattice introduces a scale (lattice unit a), 
classical scale invariance is broken and no scale invariant (instanton) solutions are 
expected to exist. Of coarse, if the radius p of the instanton is much larger than 
the lattice unit a, the values of the continuum solution in the lattice points define 
a configuration which is 'almost' a solution of the lattice Euler-Lagrange equations. 
The action of this quasi-solution will be 

A = 4n\Q\(l + c^ + 0(^)). (144) 

The size and sign of c depends on the details of the form of the lattice action. 
Depending on the sign of c, large instantons shrink or grow to get a smaller action. 

An even more serious problem is the possible existence of dislocations, topologi- 
cal artifacts. These are small, 0(a) objects which, especially when using the geomet- 
rical definition of the topological charge, might contribute to the topological charge 
and have an action below 47r|Q|. The danger is that these low-action, non-physical 
objects dominate certain topology related expectation values (like < Q 2 > /V, 
where V is the volume of the system) and might even spoil the continuum limit 
completely |2iJ. 

There are many high energy physicists who believe that fluctuating instanton 
configurations play an important role in QCD, especially for questions where chiral 
symmetry is relevant. If a discretized action describes the classical topological prop- 
erties of the theory poorly, then it might happen that continuum quantum physics 
will be reproduced on very fine lattices only. This would increase the numerical 
difficulties which are not small anyhow. 

The FP actions offer a solution to this problem [g |l|, |7| ||, g (H). We 
show now that the Euler-Lagrange equations of the FP action have scale invariant 
instanton solutions. We demonstrate first the following statement 0: 

If the configuration R on the coarse lattice satisfies the Euler-Lagrange equations 
of the FP action and it is a local minimum of A FP (R), then the configuration 
S m i„(R) on the finer lattice which minimizes the r.h.s. of the FP equation eq. ( |90| ) 
satisfies the Euler-Lagrange equations of motion as well. In addition, the value of 
the action remains unchanged, A FP (S min ) = *4 FP (R). 

This statement is easy to show. Since R is a solution of the classical equations 
of motion SA FP /<5R = 0, then the configuration S = S m j„(R) should satisfy the 
equation 

for a ny r ig. Here X nB is the sum of the four S„ spins in the block, eq. (Q). 
Eq. ( |145|) follows by taking the variation S/SH of the FP equation, eq. (|)0j). Since 
we assumed that the configuration R is a local minimum of A FP , we have A„ B > 0. 
On the other hand, the ~ n a term in eq. ( |90| ) 

- Kff ^(R„ B S„ B -|E„ B |), (146) 

takes its absolute minimum (zero) on the configura tion S m ;„ (R) satisfying eq. ( |145| ) . 
Since S TO j n (R) is the minimum of the sum of eq. (|l46| ) and „4 FP (S), it follows that 
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j(R) is a stationary point of A FP : 
^ SA FP (S) 



ss 



= 0. 



(147) 



\S=S r . 



Since eq. (146) is zero if eq. (145) is satisfied, we get A FP (R) = A (S m i n ), as we 
wanted to show. 



According to this statement, if A FP has an instanton solution of size p, then 
there exist instanton solutions of size 2p, 2 2 p, . . . with the value of the action being 
exactly 47T for all these instantons (scale invariance) . The value 4-7T follows from the 
fact that very large instantons are smooth on the lattice and then any valid lattice 
action gives the continuum value. 

It is important to observe that the reverse of the statement is not true: if the 
configuration S is a solution, then the configuration R, where R„ B = S nB /|S nB |, 
(S„ B is defined in eq. (Eq)) is not necessarily a solution. The proof fails because 
for this configuration R the minimizing configuration is not necessarily equal to 
S itself. Really, consider eq. ( |90| ) with this R configuration. We have to find 
S m in minimizing the r.h.s.. It is easy to see that S itself (which served to produce 
the configuration R) is a minimum. It is, however, not necessarily the absolute 
minimum. One can study the problem in every detail numerically |37| |39| . It turns 
out that if the configuration S has an instanton with p ^> a, then for the blocked R 
in eq. (^) S itself is the absolute minimum. In this case, R is also a solution of the 
Euler-Lagrange equations of motion and describes an instanton with radius p/2. If, 
however, p/2 is smaller than a (actually, smaller than ~ 0.7a for the FP defined by 
the block transformation eq. (|90|)) a competing, distinctly different configuration 
becomes the absolute minimum. The configuration R will not be a solution of 
the classical equations of motion, it will not describe an instanton anymore. The 
instantons with radius smaller than ~ 0.7a 'fall through the lattice', they do not 
exist. The value of the action is exactly 2tt until p is larger than ~ 0.7a. At the 
moment the instanton 'falls through the lattice', a new minimum takes over, a non- 
analyticity (a knick) occurs, and the value of the action begins to decrease. At the 
same moment, the topological charge drops to zero. This non-analyticity has no 
direct relation to the locality properties of the action. The action is local. Fig. [| is 
the result of a detailed numerical study in the 0(3) non-linear a- model |3^J. The 
existence of scale invariant instanton solutions, and the process of a small instanton 
'falling through' can be studied analytically on the example of the quantum rotor 
|4(i[ ]. The basic variables are angles in this case and the non-analyticity discussed 
above comes through a mod(27r) prescription in the FP action. The FP action is 
local, in this case even ultralocal (nearest neighbour) pp[ . 

We remark that in SU (2) Yang-Mills theory the smallest instanton the lattice 
can support has a size of ~ 0.7a quite similarly to the 0(3) tr-model case. 

To prove the statement above we only had to show that the blocking kernel 
ri°°^(R, S) in the FP equation, eq. (|S9|), is zero at its absolute minimum in R. 
This is true in all the cases we considered (cr-model, Yang-Mills theory) simply by 
construction: the normalization condition of the kernel reads 

J L»Rexp(-/3T (T (R, S)) = 1 , (148) 

which implies in the (3 — > oo limit that minR. T^^ (R, S) = 0. This can be seen 
explicitly in eq. (|9(]) and eq. (^lj) for two different RGTs in the a-model, and in 
cqs. (H,||) in the Yang-Mills theory. 
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Figure 4: Actions and charge of instantons with radii of the order of one lattice 
spacing 

The way this theorem works for the instanton solutions in the 0(3) a-model and 
in the SU(2) Yang-Mills theory has been tested in detail numerically. 



9.2 Fixed point operators 

Until now we have discussed the construction of the FP action only. We can add to 
the action different operators O; multiplied by infinitesimal sources Cj , and consider 
the RGT 

exp[-/3U'(R)-ciOi(R)-4^(R)-...] = (149) 
J DSexp[-^(A(S)+T CT (R,S)) -ci Oi(S)-c 2 2 (S)-...] . 



Start with the FP action ^4 FP (S) and take the — > oo limit. Using eq. (89) we 
obtain 

c[ 0[ (R) = c 1 1 (S mm (R)), c 2 2 (R) = c 2 2 (S roin (R)), .... 

(150) 

The FP form of the operator is reproduced by the RGT: 

Ai Of P (R) = Of p (S mm (R)) , A 2 Of P (R) = Of p (S mm (R)) 

(151) 

with Ai = c^/cj. The analysis above is similar to the general discussion in Sect. 3.4, 
except that here we also allow operators which one can not consider as part of the 
action (for example local operators which are not summed over the lattice points). 
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The eigenvalues A are determined by simple dimensional analysis if the FP is the 
Gaussian FP (which is the case for free fields and for AF theories), Sect. 8. 



Let us consider examples. Denote the field variables of a free scalar theory in d 
dimension by cf> n and Xn B 011 t ne fine and coarse lattice, respectively. Consider a 
RGT with scale factor 2. The FP field at the point n B is constructed out of the x 
fields in the local neighbourhood of n B and has the transformation law 



d-2 



2-— <p* r (x;n B ) = <p* r {4> m in(x);n = 2n B ), (152) 

where <p m in (x) is the minimizing configuration of the quadratic FP equation of the 
free field theory. 

As a second example, consider a current Jfj,(n) in QCD. The FP current is a 
local combination of the fermion fields connected by the product of gauge matrices 
to assure gauge invariance. Beyond basic symmetries (which depend on the type of 
current we consider) the details are determined by the equation 

2_3 ^(XX;^ / ;™-b) = J»{i>st{x, U m in),4> s t{x, U mm ), U min (V);n = 2n B ) , 

(153) 

where U m in is the minimizing configuration of the FP equation, eq. (|9^), while 
'ipstjipst are the configurations which make the r.h.s. of eq. (|100[) stationary. 



Only a few cases were treated explicitly until now: the FP free field eq. (152) 
and the FP Polyakov loop further the FP topological charge J37, 38, 40| |, which 
will be discussed in the next subsection. 



9.3 The FP topological charge 



As we discussed before, there are many ways to discretize the continuum definition 



of the topological charge, eq. (141). One might replace the derivatives by finite 
differences taking care of the basic symmetries of the expression. This is the so 
called 'field theoretical' definition |34j. The corresponding topological charge on 
the lattice will not be an integer number. This creates immediately a complicated 
renormalization problem since there will be contributions to this operator even in 
perturbation theory |35) , An alternative possibility is to construct an operator on 
the lattice which gives an integer number on any configuration ('geometric defini- 
tion') (36) and so it is protected in perturbation theory. This definition assigns, 
however, a non-zero topological charge to different topological dislocations of size 
0{a) which have nothing to do with instantons and which heavily distort the results 
for quantities related to topology |p9|| . 

The FP topological charge used together with the FP action avoids these diffi- 
culties. Since Q is a dimensionless quantity it satisfies the FP equation: 

Q FP (IL) — Q FP (S mm (R)) . (154) 

One can solve this equation the following way. Let us denote the minimizing con- 
figuration associated with R by s|^] n (R) (= S m i„(R) in eq. (154)). Consider now 



min 

configuration S^] n = S^„(S^-„(R)) lives on a lattice whose lattice unit is a factor 



the FP equation, eq. (|8S|), with S min on the l.h.s . The corresponding minimizing 
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of 2 2 smaller than that of the lattice of R. Iterating eq. (154) this way, we get 




Q FP (R)^Q FP (sli\ n (K)), (155) 

where S^ n (R) lives on a very fine lattice (containing big instantons with small 
fluctuations) if j is large. On such configurations any lattice definition Qi a tt gives 
the correct result 

Q FP (R)= lim Q,.«(S^(R)). (156) 

In practice, the j ; = 2 value agrees with the j = oo limit if the geometrical definition 
is used |37|, p9[ fo r Qiatt- The construction of goes the same way in Yang-Mills 
theories. 



9.4 There are no topological artifacts if FP operators are 
used 



We show now that using the FP action and the FP topological charge, the action 
of any configuration (classical solution or not) satisfies 

A fp (R)>4tt\Q\, (157) 

like in the continuum classical theory. There are no configurations with topological 
charge whose action is below that of the instanton solutions. 

In the equation A FP (R) = A FP (S m i n )+Ta\'R.,S m i n ) we can express A FP (S m i n ) 
again in terms of the sum of the action and the kernel on the next finer lattice, 
etc. Iterating further we get arbitrarily close to the continuum, where the inequal- 



ity eq. (157 ) is certainly satisfied. Since the kernels are non-negative (Sect. 9.1), 



cq. (157) follows. 



10 Chiral symmetry on the lattice 



Although it has not been demonstrated rigorously, it is strongly believed that chiral 
symmetry is spontaneously broken in QCD. Without this dynamical effect QCD 
can not describe the observed properties of low energy hadrons. The formal order 
parameter of chiral symmetry (^){x)tp{x)) is equal to the trace of the quark propa- 
gator averaged over the gauge field configurations. The trace is proportional to the 
quark mass m q , hence in the m q — > chiral limit we can get a non-zero value only 
if there is an appropriate accumulation of small eigenvalues of the massless Dirac 
operator p2| . Actually, there exists a rigorous theorem |Q concerning the zero 
eigenvalues of the eigenvalue equation in the continuum (we consider one flavour, 
N f = l) 

1 ^-iA tJ ,{x))^ x) {x)=\^ x \x). (158) 

According to this theorem, if the background gauge field has a non-zero topological 
charge Q, then there necessarily exist zero eigenvalues. The corresponding eigen- 
vectors are 75 eigenstates: 75^^ l = l- ^ we denote by tir and the number 
of right and left handed eigenvalues, the theorem says 

n-L - n R = Q . (159) 
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This result might have important consequences on the low energy properties of 
QCD. A possible intuitive picture of a typical gauge configuration in QCD is that 
of a gas or liquid of instantons and anti-instantons with quantum fluctuations p4| . 
For such configurations, the Dirac operator is expected to have a large number 
of quasi-zero modes which could be responsible for spontaneous chiral symmetry- 
breaking. 

Typical gauge field configurations contributing to the path integral are not 
smooth, however. In addition, the regularization which is necessary to define the 
quantum field theory breaks some other conditions of the index theorem as well. 
Standard lattice formulations, for example, violate the conditions of the index theo- 
rem in all possible ways: the topological charge of coarse gauge configurations is not 
properly defined and the Dirac operator breaks chiral symmetry in an essential way. 
Close to the continuum some trace of the index theorem can be identified |f5|, [46| . 
By improving the chiral behaviour of the Wilson action a strict connection between 
the small real eigenvalues of the Dirac operator and the geometric definition of the 
topological charge can be established even on coarse configurations j47fl . In general, 
however, the expected chiral zero modes are washed away and, even worse, the real 
modes which occur create serious difficulties in simulations. They lead to 'excep- 
tional configurations' where the quark propagator becomes singular when the bare 
mass is still distinctly different from its critical value Q] . 

Even in the continuum limit, the explicit chiral symmetry breaking in the action 
leaves its trace behind in the quantum theory creating technical difficulties. The 
quark mass receives an additive renormalization, hence the bare quark mass should 
be tuned to its (/? dependent) critical value to get the pion massless. Against 
the expectation that a partially conserved current is not renormalized, the axial 
vector current suffers also from renormalization. Further, there is mixing between 
operators in nominally different chiral representations creating a special difficulty 
in weak matrixclcmcnt calculations. 

Fixed point actions offer a bold solution to all these problems. The good chiral 
behaviour of the FP action follows essentially from the fact that its fermion part 
satisfies the Ginsparg- Wilson remnant chiral symmetry condition |j49fl , for which 
no solution was known previously. This is the point, where the FP action meets 
a seemingly unrelated approach to overcome the problems of lattice regularization 
concerning chiral symmetry |30[ ^l) . The purpose of this section is to discuss these 
recent developments. 



10.1 'On-shell' symmetries 

As the example of instanton solutions illustrates, scale invariance, which is obvi- 
ously broken by the lattice, remains a symmetry on the classical solutions of the 
FP action. This observation can be generalized [j52| . For any symmetry of the 
classical countinuum theory (arbitrary translation, rotation, chiral rotation, etc.) a 
representation can be defined on the classical solutions. 

For the arguments, let us use the language of the non-linear cr-model. Assume 
that the confuguration R, which lives on a lattice with lattice unit clq , is a solution 
of the Euler-Lagrange equations of motion of the FP action. As we have discussed 
in Sect. 9.1, the minimizing configuration S^- n (R) associated with R is also a 
solution. This configuration is defined on a lattice with unit a\ — o.q/2. We can 
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iterate this process: to the solution S^' n (R) there corresponds again a minimizing 
configuration S^ n (S^ n (R)) which lives on a lattice with lattice unit a 2 = ao/2 2 , 
etc. Iterating a large number of times we get arbitrarily close to the continuum 

lim S«„ = S(z,R), (160) 

j— »oo 

where S(x, R) is a solution in the continuum. On this solution all the symmetry 
transformations of the continuum theory are well defined. Take, for example an 
infinitesimal translation: S — > S e where S e (x,R) = S(x + e,R). The configuration 
S e is also a solution which we can block back to the starting lattice with unit 00 
producing a a new solution R e there [J The process 

R„ -> S(x,R) -> S(x + e,R) -> R^ , (161) 

defines a representation of the group of continuous translations. It goes similarly 
for other symmetry transformations of the classical continuum theory. 

We emphasize that the procedure above makes sense only for configurations 
which are solutions of the classical equations of motion. Therefore the name: 'on- 
shell' symmetry. 



10.2 The index theorem on the lattice (I) 

Consider the eigenvalue equation of the FP Dirac operator 

Y,h FP {U) n ^) . (162) 

n' 

The index theorem refers to the A = eigenvalues, in which case the eigenvalue 
equation above is reduced to classical Euler-Lagrange equation (Dirac equation). 
Due to the 'on-shell' symmetries discussed in Sect. 10.1, we suspect that the FP 
action will lead to an exact index theorem on the lattice independently whether the 
background gauge configuration is smooth, or coarse. This is really so |5^]. The 
proof is based on the FP equations for the gauge action (eq. (j92|)) and for the Dirac 
operator (eq. ( |100| ). 

Consider a gauge configuration V on a lattice with lattice unit a>o- Assume that 
V has a topological charge Q. As discussed in Sect. 9.3, the FP topological charge is 
defined through a sequence of minimizing configurations U^} n (V), U^l n (U^} n (V)), . . 
Here Uzj n lives on a lattice with lattice unit aj — a /2 J . By construction all these 
configurations have the same topological charge Q, which is the charge of U^°l in 
the continuum, at the end of the sequence. 

Assume that the FP Dirac operator h FP has a zero eigenvalue for the background 
field V 

J2 hFP ^ B ,n' B XW B =0. (163) 

n B 

6 As discussed in Sect. 9.1 blocking a solution might not lead to a solution again if we land 
on a very coarse lattice. Assume, the original solution R is an instanton with radius p slightly 
larger than the critical size ~ 0.7a. Blocking back the corresponding continuum solution after 
translation, we get a translated instanton with the same p on the coarse lattice. Since the critical 
size depends slightly on the position of the instanton on the lattice j37| ^5sj] , the translated small 
instanton might 'fall through the lattice'. To avoid this complication, assume that the solution R 
is away from this critical region. 
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Eq. ( |163[ ) is the Dirac equation, i.e. the classical Euler-Lagrange equation of motion. 
We proceed now as we did in the case of the classical solutions in Sect. 9.1. Taking 
the variation of eq. (100) with respect to \n B and using eqs. ( |96] , |163| ,86) we find 
that ipn' which makes the r.h.s. of eq. (100) stationary, satisfies the Dirac equation 
for the background field uj^J (V). Since eq. (10C) is quadratic in the fermion field, 
the stationary configuration is unique. Hence the opposite statement is also true: 
if ipn is a solution of the Dirac equation for the background field U^J n (V) then 
Xn B = bfJ2n LJ (^)n B ,n' l Pn is a solution on the coarse lattice for the configuration 
V. 



This result can be iterated towards finer lattices. Having a zero mode eigenfunc- 
tion ip!^ for U^J n there is a zero mode on the next finer lattice for the configuration 

(2) 

Umin> e ^ c - At every step we can also move in the opposite direction. Going up, 
we finally get to a very smooth configuration for which the Atiyah-Singer theorem 
is applicable. This way we get that for each gauge configuration of the sequence, 
the FP Dirac eigenvalue problem has the same number of zero modes and with the 
same helicities as in the continuum. The statement on the helicities follows from 
the relation between the eigenvectors of the z-th and (i — l)-th step 

C = h E^™L)™<^>^<.) . ( 164 ) 

and from the fact that lj is trivial in Dirac space (Sect. 5.3). 



We have to emphasize that for the validity of the index theorem we need a 
classically perfect regularization both in the fermion and in the gauge sector. This 
is also obvious from the fact that the index theorem connects a quantity defined in 
the gauge sector (topological charge Q of the gauge configuration) with that coming 
from the fermion sector (the number of R and L zero eigenmodes) . 



10.3 Free FP fermions 



A good starting point to understand the chiral properties of the FP action is the 
case of free FP fermions. The FP of free theories (scalar, gauge, or fermions) 
can be constructed analytically and their properties can be studied in every detail 
fill ||li, Sect. 8. 

Write Urnin = 1 in the recursion relation eq. (Bq) and consider the example of 
the simple block transformation ui(2tib — n) = 2~ if n is in the hypercube whose 
center is indexed by and zero otherwise, Sect. 5.3 . Starting with the Wilson 
action (with m q = 0) on a very fine lattice, the iterated RGT drives the system to 
a FP. In Fourier space the FP propagator reads M, [l2| 

Here q is defined over the Brillouin zone. The poles of this propagator give the exact 
continuum spectrum, as expected. There is no doubling. (The massive doublers run 
out to infinity under the repeated RGT.) The action is a periodic analytic function 
of q, in configuration space it is local. Complying with the no-go theorem |p4f , the 
propagator, and so the action, break chiral symmetry. This symmetry breaking is, 
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however, very special: 

{h- 1 ^ 5 } = — 7 5 . (166) 

The chiral symmetry breaking term 2/ Kf in the propagator comes entirely from the 
block transformation. The effect of the chiral symmetry breaking in the propagator 
shows up only at n = in configuration space. No physical properties are influenced 
by that |. 



10.4 The Ginsparg- Wilson remnant chiral symmetry condi- 
tion 



More than fifteen y ears ago, Ginsparg and Wilson |49| suggested to elevate the 
free field result eq. (|l66| ) into a general principle which defines (remnant) chiral 
symmetry on the lattice for interacting theories also. The condition is derived in 
Ref. pfij] by considering RG transformations starting out of the continuum theory. 
Since, however, in the derivation the assumption is made that the blocked action is 
quadratic in the fermion fields (which is true only in a free field theory) it might be 
better to consider the GW remnant chiral symmetry condition as an genial, intuitive 
suggestion: the lattice regularization is expected to preserve the consequences of 
chiral symmetry in interacting theories also (in particular in QCD) if the Dirac 
operator satisfies 

\{Kl(U),f} = R m AU)i\ (167) 



or equivalently 



\{h nn > ,7 5 } = i^ 5 •Rh) nn , . (168) 



Here /i([/) _1 is the quark propagator over the background gauge field U and the 
matrix R nn > is trivial in Dirac space. Both h(U) nn > and R(U) nn i are local. The 
locality of R in the remnant chiral symmetry condition is absolutely essential. It is 
also highly non-trivial, since ft -1 is a non-local quantity. 



Ginsparg and Wilson have shown |49[ that eq. (167) implies the correct triangle 
anomaly on the lattice. In the presence of gauge fields, however, no solution was 
found by the authors and the paper remained practically unnoticed for a long time. 



10.5 Solutions of the remnant chiral symmetry condition 

It has been noticed recently that the FP action of QCD satisfies the remnant chiral 



symmetry condition eq. (167) [|55| |56|]. This result is a direct consequence of the FP 
equations eqs. ( |92| , |98| ). 

Eq. ( p8| ) can be solved by iteration. If the scale change s of the transformation 
is (infinitely) large, the fixed point can be reached in a single renormalization group 
step. In this case, the fine field U m i n lives on an infinitely fine lattice, and it is 

7 At Kf = oo the FP action becomes chiral symmetric, but at the same time, non-local [^| in 
consistency with the no-go theorem |24[ 
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smooth. The Dirac operator h over this field is arbitrarily close to the massless 
continuum Dirac operator. Since to is trivial in Dirac space, the second term on the 
r.h.s. of eq. (pq) is chiral invariant, which leads to the remnant chiral symmetry 



condition in eq. (167) with R nn > = ■ S nn i. If s is finite, one can start again on 
an infinitely fine lattice, but the transformation should be iterated. One obtains in 
every iteration step some contribution to the chiral symmetry breaking part of h , 
and a non-trivial, but local R builds up in this case. 



Very recently, another solution of eq. (167) was presented |57|| . This solution 
grew out the overlap formalism |5l], |5^] and seems to be unrelated to renormalization 
group considerations. The solution has the structure h(U) n _ n i = 5 nn i + V(U) nn > 
(only the space indices are indicated explicitly), where the unitary matrix V has 
the form 

V = X^J- . (169) 

Here X(U) n<n i, X* = 7 s A7 5 is an ultralocal interaction, whose form is explicitly 
given |57| . In spite of the square-root in the definition of V, the action is obviously 
local in the free case (U — 1), and is, presumably, also local in general. It is 
a little exercise to show that the Dirac operator h defined this way satisfies the 



remnant chiral symmetry condition cq. (167) with R n .n' = l/2S n .n'- An advantage 



of this solution is that it is explicitly given, while to find the FP action the defining 
classical equations have to be solved. A disadvantage is that it is significantly 
less shortranged than the optimized FP Dirac operators. They decay exponentially 
with a constant 7 = 0(3) in the exponent, while the solution discussed above decays 
only with 7 = 0(1). Another difference is that the overlap formalism addresses the 
issue of chiral symmetry only, while the FP action is classically perfect. For free 
fermions, for example, the action defined by eq. ( |165| ) has an exact spectrum, while 
the spectrum of the Dirac operator above is similar to that of the Wilson action. 
Nevertheless, it is interesting that the structure of eq. ( |169| ) alone assures the GW 
condition which might be a very useful observation in itself. 



10.6 Currents 



In the following sections we shall need different currents which we construct here 
and study some of their properties. Consider a fermion action which is quadratic 
in the fermion fields 

A f $ t i/),U) = i>m,<h ml ^{U)xlj n , . (170) 

m' ,n r 

Most of our discussion will be general, in particular we do not assume that Af is a 
FP action, or that it satisfies the Ginsparg- Wilson relation. 

A) U(l) gauge field, vector current, sum rules 

We shall use the example of the electromagnetic current to illustrate the pro- 
cedure. Consider a U(l) gauge field (electromagnetic field) coupled to a charged 
fermion and construct the current which belongs to the global vector symmetry. 
Following the usual steps of defining a Noether current we consider the local trans- 
formation 

Sip n = ie n ip n , Sip n = -ip n ie n , (171) 
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and obtain 



SAf = -i^ l € n 7Mn), (172) 

n 

where 

li .III 4>rn> - l/j m >h(U) m ',n1pn] . (173) 



The derivative is defined as 



Vft/Cn) = f(n) - f(n - fi) . (174) 



We want to find a vector field Jfi(n) which satisfies eq. (173). A solution is 



Jfi{n) = "I V i>m< , , , ■, h(U) m ',n'1pn' , (175) 

where A^(n) is the vector potential, U^{n) = e lA ^ n \ The action .4^ is gauge 
invariant: ip m i and ?/v are connected by a product of U fields along paths running 
between the fermion offset m',n'. Consider a fix fermion offset m',n'. Since 

U u (m) = iS^vSnmU^ri) (176) 



SA M (n 

we get that every path (between m',n') going through the link (n,n + fi) gives a 
contribution to J^(n) which is equal to the contribution of this path to h m ^ n > itself, 
times (s + — s_), where s + (s_) is the number this path runs through (n, n + fx) in 
positive (negative) direction. This U(l) current is obviously gauge invariant. 



We want to show that the current in eq. (175) satisfies eq. (173). Consider 

V M J M (n) = (Ji(n) - Ji(n - 1)) + (J 2 (») - J 2 (n - 2)) + . . . . (177) 

To Ji (n — 1 ) all the paths contribute which go through the link (n — 1 , n) . Take first 
n m' , n n' . All the paths which run into the point through the link (n — 1, n), 
(and so contributing to J\(n — 1)), will run out of this point through one of the 
links (n,n + v) (and so contributing to J^n)), etc. These contributions cancel. A 
similar consideration shows that if n — to', or n = n' the contribution to V M J^(n) 



is h n , m i and — /i m ',n, respectively, leading to eq. (173) 



If the action .A* is invariant under the reflection of one of the coordinate axis, 



then J^{n) in eq. (175) will transform like the vector potential A„(n) under this 



transformation, i.e. like a lattice vector field, as it should. 



Consider the current with U = 1. The current is translation invariant in this 
case and we can write 

J m(") = ^2 $™' T n( m ' -n,n' - n)ip n , . (178) 

m 1 ,n' 

Although the form of the current depends on the details of the Dirac operator, 
the following sum rules, which test the long-distance properties of the current, are 
generally valid if the action preserves cubic symmetry: 

5^r p (m,m-0 = -i P Ai, (179) 
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£ m„r M (m, m - I) = JJk+M hl , (i 80 ) 

m 

where we introduced the notation h(U = l) m „ = h m - n . Since these sum rules 
enter the derivation of the axial anomaly ^, we sketch their derivation. 

Fix the fermion offset m! — n' — —I. Using translation symmetry, we can take 
m! = 0, n' = I. Summing over n 

^)r M (-n,i-n), (181) 

n 

we have to consider all the paths in h^i and for each path count the links going 
in the ±/t direction with the weight ±1. Since the paths run from (0, 0, ... , 0) to 
(h, I2, ■ ■ ■ , Id), each path gives leading to 

J2~Tn(-n,l-n) = l fi h- l (182) 

n 

Taking I —I and changing the summation variable — n — > m we get the first sum 
rule in eq. (179). 

Consider now 

J2n v T„(-n,l-n). (183) 

n 

Due to hypercubic symmetry, a path and its reflection to the point in the middle 
of I contribute the same amount to h-i. Take the contribution of an arbitrary path 
and of its central reflected path to eq. ( |183| ). Consider first (j, ^ v. A path going 
through (n, n + /t) contributes n v . The reflected path gives l v — n v , so the average 
is 1 V I1. Summing over all the links in the /t direction gives l^l v /2. Consider now 
\i = v. The link (n, n + (1) on the path contributes in^ depending on the direction 
of running through. Summing over all the links gives 1^(1^ — l)/2. We obtain 

nuT^-n, l-n)= W • (184) 

n 

Taking I — > — I and changing the summation variable —n — » m we get the second 
sum rule in eq. (180). 

B) QCD flavour currents 



Consider eq. (170) with Nf flavours and assume that h is flavour symmetric. 
Denote the flavour generators by r a , a = 1, . . . , N'j — 1 with 

[r a ,r b } = if abc r c , tr(r a T b ) - l -8 ab . (185) 

The fermion action has a U(Nf) flavour vector symmetry. As we discussed before, 
the axial symmetry is broken, even if we set the bare mass to zero. As usual, the 
axial currents are constructed from the chiral symmetric part of the action using 
^■sym = 5[^:7 5 ]7 5 - There is a certain freedom in constructing the vector currents. 

8 In Rcf. J49| the current was constructed in a particular way, and the sum rules above were 
shown explicitly to be satisfied. The particular form of the current in Ref. |flSf| is not convenient 
in different theoretical considerations, however. 
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Using the full action, conserved vector currents are obtained. On the other hand, the 
asymmetry between the vector and axial case generates strange terms in current 
algebra relations and it requires extra work to show that they go away. Let us 
construct both currents with /isym- Using the steps we discussed in A) we get 

V^(n) = ^ n T a (h SYM ^) n - (# SYM )„r>„ , (186) 

V„A«(n) = -^ n r a (~/ 5 h SYM iP) n - (#sym7 5 )„ r a ^ n . (187) 
None of these currents are conserved if the equations of motion hip = is used. 



In order to find the currents themselves we introduce the flavour gauge matrices 
W fl (n) = l+iw^ (n)T a + . . . and extend the product of colour U matrices in eq. (170) 
along the paths between the fermion offsets m',n' by the product of W matrices 
along the same paths. Then the vector current is defined by 



TO 



8w a (n 



■(hsYM)m'n'{U,W)lpr, 



(188) 



This current satisfies eq. (186) as can be shown using the same arguments we used 
in A). Eq. (188) can also be written as 



(189) 



where T ^ is trivial in flavour space. Similar considerations lead to the axial current 



^» = E ^m'r a Tl(m',n';n;U)^ 



(190) 



where 



r 5 — r -y 5 — -T 5 r 



(191) 



The flavour singlet vector and axial currents are obtained from eqs. (186 187) by 
replacing r a by the unit matrix in flavour space. 



10.7 The index theorem on the lattice (II) 

In Sect. 10.2 we demonstrated that the index theorem is valid on the lattice if the 
FP action is used. The method was to connect the lattice problem with that in 
the continuum (where the index theorem is true) using the FP equations. One can 
prove, however, the index theorem directly on the lattice without refering to results 
obtained earlier in the continuum |5(| . We shall discuss this derivation briefly since 
it leads to additional results and insight. 

We shall suppress the index FP on the fixed point Dirac operator h FP in the 
following. To avoid some singularities at intermediate steps of the calculation we 
introduce a small quark mass which will be sent to zero at the end: 

(U) —> h m 'n' (U) = h m i n > (U) + m q 5 

m'n' ■ 

(192) 

Due to eq. ( |168| ) the chiral symmetry breaking part of h satisfies the relation 

-{hm'n'-, 75 }75 = {ht Rh) m > n > +m q 8 rn , n , . (193) 
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This equation plays a basic role in deriving the results. 



Consider the flavour singlet axial vector current (take Nf — 1) and write 

ftsYM = ^[ft,75]75 = h- ^{h, 75)75 . (194) 



Using eq. ( |193| ) we obtain 

-Va.75 {h ] Rhi£) n + (^Rh) nl5 i, n (195) 
-2m q -0„7 5 V'„ ■ 

The first term on the r.h.s. is proportional to the equation of motion hip = 0, 
the second comes from the remnant chiral symmetry condition, while the last term 
is due to the explicit symmetry breaking. Consider the expectation value of this 
equation over a background gauge field. For any operator O this expectation value 
is defined as 



1 

Z~U) 



(0(iP, ip, U)) = / DiPDipO(iP, V, U) exp{- ^ ip m h mn (U)ip n } , 



(196) 



where Z(U) is given by the analogous integral without the 0(ip, ip, U) factor. The 
contribution from the first term is zero. Indeed, 

(197) 



{ipnlb{hip) n ) = -777777 / DipDip i/)„7 5 - rr - exp{- ) ipm'h m > n >ip n >} 
Z(U) J 5lp n ^ 



in .11 



and partial integration by ip gives then tr7s = 0. Integrating out the fermions in the 
expectation value of the third term in eq. (195) gives 2m q ti{^h~^ l ) , where the trace 
is over colour and Dirac space. We insert here the closure relation of eigenvectors 
of the hermitian matrix H = /175. It is easy to show using the Ginspar- Wilson 
relation that the zero eigenvalue eigenfunctions of H (if there are any) are chiral 
and they are zero eigenvalue eigenfunctions of h as well. In the m q — > limit only 
these eigenfunctions contribute leading to 

"L "R 

lim m q (^„7 5 -0„) =^^*(n)<ip£\n) - V ip^* {n)ip^ (n) , 
mq 3 =i i=i (198) 

where ip S (n) , tp^ (n) are the normalized A = right- and left-handed eigenfunc- 
tions of the eigenvalue problem of the Dirac operator. Finally, integrating out the 
fermions in the expectation value of the second term in eq. (195) and summing over 
n we obtain 

= 2 TriTshR) + 2{n R - n L ) , (199) 

where the trace Tr is over colour, Dirac and configuration space. The first term is a 
pseudoscalar gauge invariant functional of the background gauge field. For smooth 
fields U its density is expected to be proportional to e MI/a/3 F^ v (n)F®p(n). Really, 
one can show |56| 

- $ n7s (h^Rhip) n + {^Rh) nl5 ip n ) - ^^F;(n)^(n) 

dZ7T (200) 
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up to terms which are of higher order in the vector potential and/or derivatives. 
This result follows from Ref. where the anomaly of the flavour singlet axial 
current was derived using the remnant chiral symmetry condition and the sum 
rules on the currents we discussed in Sect. 10.6. Further, it can be shown that the 
first term in eq. (199) is a FP operator: it is reproduced under a RGT. It follows 
than that 



Tr( 75 W?) = Q 



FP 



(201) 



where Q is the FP topological charge defined in the Yang- Mills theory, Sect. 9.3. 
Eq. (201) is a remarkable connection between the gauge and fermion parts of the 
FP theory. 



10.8 The spectrum of the Dirac operator satisfying the rem- 
nant chiral symmetry condition 



Let ip be a normalized solution of the eigenvalue equation, hip = Xip. Eq. (168) 
(written in the form ^(h + /i^) = h^Rh) implies ReA = |A| 2 (-0, Rip). For the special 
case R nn ' = • 8 nn i the eigenvalues A lie on a circle of radius /tf/2 touching 
the imaginary axis, as shown in fig. |^. In general, the eigenvalues lie between two 
circles tangent to the imaginary axis - the circle described previously and a smaller 
one as shown in fig. ^. The smaller circle has a finite radius which follows from 
the locality of R. This property excludes exceptional configurations, appearing e.g. 
for the Wilson action and causing serious problems in numerical simulations. As a 
consequence of h) = hj^h the solutions with ImA ^ come in complex conjugate 



pairs and satisfy (V'jTsVO = 0. Due to eqs. (167,168), those with A = are chiral 
eigenstates, 75^ = ±ip- 




Figure 5: The spectrum of the fixed point Dirac operator: a) for the special case 
Run' = l//tf • <5nn' the eigenvalues lie on a circle, b) for a more general set of block 
transformations the eigenvalues lie between the two circles. 



10.9 Lattice QCD without tuning, mixing and current renor- 
malization 



As we discussed at length before, lattice regularized, local QCD actions break chiral 
symmetry under general conditions. Wilson fermions have a dimension five symme- 
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try breaking operator whose effect on the physical predictions goes to zero in the 
continuum limit. It leaves, however its trace behind in the form of additive quark 
mass renormalization, axial current renormalization and mixing between operators 
in nominally different chiral representations. There is a significant recent progress 
in calculating these renormalizations in a theoretically controlled non-perturbative 
way [[59]]. Nevertheless, the situation is not really pleasing theoretically and the 
technical difficulties are also significant. 

On the other hand, if the Dirac operator satisfies the remnant chiral symmetry 
condition, then all the physical consequences of chiral symmetry survive the lat- 
tice regularization without quark mass tuning, without the mixing of operators in 
different chiral representation and without flavour non-singlet axial current renor- 
malization We do not discuss here the details of the arguments (which are 
based on chiral Ward identities following the seminal paper Ref. ]6lj| but mention 
the intuitive reason only. Consider the divergence of the flavour non-singlet axial 
vector current for m q — 

V M v4;» = ^„T a 7 5 (/«/>)„ + 0M)„75T Q V>n 

~^ n T a (hrftRhil>) n + {4>h lb Rh) n r a i> n . (202) 

The first term on the l.h.s. is proportional to the equations of motion. Observe 
the special structure of the chiral symmetry breaking second term: it contains two 
h factors. If the divergence of the axial vector current, and so this term enters a 
Green's function these two h factors will cancel the two propagators which arrive 



when the -011 and ipn fields in eq. ( 202 ) are paired with the fermion fields of the 
other operators in the Green's function. This breaking term, therefore is expected 
to produce contact terms only which will not influence the physical predictions. 



10.10 Exact chiral symmetry on the lattice 

Looking at the results of the previous sections, one might suspect that there exists an 
underlying hidden exact symmetry on the lattice which replaces the standard chiral 



symmetry in the formal continuum formulation. This is really so 62 . Luscher 
observed that the remnant chiral symmetry condition eq. ( |168| ) implies an exact 
symmetry of the fermion action, which may be regarded as a lattice form of an 
infinitesimal chiral rotation. Really, the transformation 

Sip n = ie a T a ^2~/ 5 (l - Rh) nin ri/J n > , 

n' 

5i> n = iJ2tM>{l-hR) n ,, nl5 e a T a , (203) 

n' 

which may be considered as a locally sm eared version of the standard 75 transfor- 



mation, is a symmetry of the action eq. ( 170 ). There is no contradiction with the 



no-go theorem |2j], since the symmetry is realized differently as has been assumed 
when proving the theorem. A large circle is closed by this observation. 

One might hope that the recent developments discussed in the last sections above 
will help to solve further longstanding difficulties of lattice regularization. 
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